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Abstract: Two-dimensional s[(ra) quantum Toda field theory on a sphere is considered. 
This theory provides an important example of conformal field theory with higher spin 
symmetry. We derive the three-point correlation functions of the exponential fields if one 
of the three fields has a special form. In this case it is possible to write down and solve 
explicitly the differential equation for the four-point correlation function if the fourth field 
is completely degenerate. We give also expressions for the three-point correlation functions 
in the cases, when they can be expressed in terms of known functions. The semiclassical and 
minisuperspace approaches in the conformal Toda field theory are studied and the results 
coming from these approaches are compared with the proposed analytical expression for 
the three-point correlation function. We show, that in the framework of semiclassical and 
minisuperspace approaches general three-point correlation function can be reduced to the 
finite-dimensional integral. 

Keywords: Integrable Field Theories, Conformal and W Symmetry. 



September 24, 2007 



Contents 

Introduction 1 

1. Toda Field Theory 3 

2. Differential equation 12 

3. Classical limit (heavy exponential fields) 23 

4. Classical limit (light exponential fields) 31 

5. Minisuperspace limit 40 
Conclusion 43 

A. The Coulomb integrals 44 

B. Simplification of the integral (4.18) 47 

C. Properties of the s[(3) Coulomb integral 48 

D. Useful formulae 50 

E. Example of application of the recursive relation (5.6) 51 



Introduction 

It is well known, that the problem of integrating over all metrics modulo diffeomorphism on 
a two-dimensional surface can be reduced to studying the quantum Liouville field theory 
[1]. First attempts to solve this theory were transformed into a beautiful and complete 
theory known as the two-dimensional conformal field theory [2]. This theory is exactly 
solvable because the algebra of generators of the conformal symmetry in two dimensions, 
which governs the theory, is infinite dimensional. It coincides with the Virasoro algebra, 
which is the central extension of the algebra of vector fields on a circle. It is well known, 
that Virasoro algebra can be obtained as a quantum Drinfeld-Sokolov reduction of the affine 
sl(2) algebra. The same construction can be generalized to the case of general affine simple 
Lie algebra g. As a result, after reduction one obtains associative algebra (W algebra), as 
an additional infinite dimensional symmetry consistent with conformal symmetry, i. e. as 
a direct extension of the Virasoro algebra [3]. Two-dimensional Toda field theory (TFT) 
associated with simple Lie algebra g generalizes Liouville field theory in a similar sense. 



- 1 - 



The algebra of the generators of the symmetry, which governs TFT dynamics, coincides 
with W algebra (associated with the corresponding Lie algebra g). 

Due to its geometric interpretation [4, 5], TFT is relevant in the investigation of the 
W strings and W gravity (see for example Refs [6, 7]). It also provides an important 
example of non-rational conformal field theory with higher spin symmetry and hence has 
its own interest. This higher spin symmetry manifests itself also in rational conformal 
field theories, which describe the critical behavior of many interesting statistical systems, 
like for example Zn Ising models (parafermionic CFT [8]), tricritical Ising and Z3 Potts 
models, Ashkin- Teller models and also in the large variety of integrable statistical systems 
studied and solved in Refs [9, 10]. The results derived in conformal Toda field theory can 
be applied to study of the short-distance asymptotics of the correlation functions in the 
massive integrable quantum field theory, which is known as affinc Toda field theory, as well 
as to calculation of the vacuum expectation values of the exponential fields in this theory 
(see for example Refs [11, 12, 13]). As conformal TFTs appear by the quantum Hamiltonian 
reduction of the WZNW models (see for example [14]), they can be also applied to study 
WZNW models with non-compact Lie algebras. 

There has been much progress in understanding Liouville field theory (sl(2) TFT) and 
hence in the conformal field theory itself in the middle of 90's. In particular, the three-point 
correlation function was found explicitly for arbitrary exponential fields [15, 16, 17, 18]. 
Known three-point correlation functions, together with the fact, that conformal blocks are 
completely determined by the conformal symmetry, solve the conformal bootstrap problem 
in Liouville field theory. 

Conformal Toda field theory is much more complicated than the Liouville field theory. 
One of the main reasons is that in TFT we need in general case more data to solve the 
conformal bootstrap problem [19]. In particular, this difficulty manifests itself in the fact 
that contrary to the Liouville field theory it is impossible to write down the differential 
equation for the four-point correlation function, which contains one completely degenerate 
and three arbitrary fields [19, 20]. In Liouville field theory it allows to write down func- 
tional relation for the general three-point correlation function, which in some domain of 
parameters has a unique solution (see for example [21]). In TFT this procedure fails (see 
section 2 for details). It means that other methods should be applied. It is interesting, 
that the difficulty of such a type appears already at the classical level, where the problem 
of finding the solution to the sl{n) classical Toda equation for n > 2 with three singular 
points (which determines so called "heavy" semiclassical limit of the three-point correlation 
function) reduces to the problem of studying Fuchsian ordinary differential equation with 
accessory parameters (see section 3) . Accessory parameters are absent in the Liouville case 
(s[(2) TFT) and this is the reason why this theory is rather simpler. 

This paper is the first of two papers, devoted to study the correlation functions in 
the sl{n) TFT, which can be found analytically (may be in terms of finite dimensional 
integrals). It is organized as follows: in section 1 we briefly remind some basic facts 
about conformal TFT and propose an analytical expression for the three-point correlation 
function of the exponential fields in the case, when parameters of one of the field take 
the special values (see Eq (1.39)). We give also several another examples of correlation 



-2- 



functions, which can be expressed in terms of known functions. In section 2 we present 
the derivation of the proposed three-point correlation function (1.39) by using the special 
properties of the operator algebra of degenerate fields. In sections 3 and 4 the semiclassical 
analysis of the theory is developed. In the section 3 we study the case, when all exponential 
fields in correlation function are "heavy" (i. e. have parameters proportional to the opposite 
coupling constant) and in the section 4 we study the case, when all exponential fields are 
"light" (i. e. have parameters proportional to the coupling constant). In section 5 we 
study the minisuperspace approach to the sl(n) TFT. We show, that in the case of light 
exponential fields, as well as in the minisuperspace limit, three-point correlation function 
can be expressed in terms of finite dimensional integrals. In both cases, semiclassical and 
minisuperspace asymptotic is in complete agreement with the proposed quantum results. 
The calculation details and useful formulae are given in the appendices. 

In the second part of this paper [22] we will give more detailed description of the cor- 
relation functions in conformal TFT, which can be expressed in terms of finite dimensional 
Coulomb integrals. 



1. Toda Field Theory 



We start by recalling some basic facts and notions. The Lagrangian of the 5l(n) conformal 
TFT has the form 



n~l 



(1.1) 



fe=i 



here ip is the two-dimensional (n — 1) component scalar field (p = {ipi . . . ipn-i), b is the 
dimensionless coupling constant, fi is the scale parameter called the cosmological constant 
and (ek,ip) denotes the scalar product, where vectors are the simple roots of the Lie 
algebra s[(n) with the matrix of the scalar products Kij = {ei,ej) (Cartan matrix) 



Kij = 



/ 2 -1 
-1 2 -1 
-1 .... 









Vo 



... -1 

-1 2 -1 
0-12/ 



(1.2) 



In the following we will use standart for the two-dimensional physics complex notations: 



— d = — 

dz' dz 



Z = Xl + 1X2-1 Z = Xl — 1X2, d ■ 

and introduce the notation for the measure 

(fz = dxidx2- 

Total normalization of the Lagrangian (1.1) is chosen in such a way, that 



(1.3) 



(1.4) 



(pi{z, z)(pj{0, 0) = —Sij log -|- . . . at z — > 0. 



(1.5) 
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In is useful to write TFT action explicitly in reference metric gab on a surface 

n— 1 



Atft= I (^9°'(a<.V,*»>) + ^ii + *'Ee'<='"-^' ) \rg<fx, (1.6) 



fe=l 



here R is the scalar curvature of the background metric.^ If the background charge Q is 
related with the parameter h as 

Q=i^+T]p (1-7) 



6, 

with p being a Weyl vector (half of the sum of all positive roots), then the theory (1.6) 
is conformaly invariant.^ Moreover it ensures higher-spin symmetry: there are n — 1 
holomorphic currents W^(2:) with the spins k = 2, 3, . . . ,n, which are expressed through 
the field if via the Miura transformation [23] 

n— 1 n 

Y[{qd+{K-i,d^)) = ^ W"-^^)(g5)^ (1.8) 

i=0 fe=0 

where 

q = h + l/h (1.9) 

and vectors are the weights of the first fundamental representation tti of the Lie algebra 
si{n) with the highest weight u\ (first fundamental weight) 

hk = u>i- ei ek-i- (1-10) 

In particular, it follows from Eq (1.8), that the currents W^{z) = 1, W-'-(2;) = and the 
current 

W^{z)=T{z) = --{d^f + {Q,d'^) 

is the stress-energy tensor of the theory, which ensures local conformal invariance of TFT. 
The currents W^{z) form closed Wn algebra, which contains as subalgebra the Virasoro 
algebra with the central charge 

c = n-l + 12Q2 = (n - 1)(1 + n{n + 1)(5 + 5"^)^). (1.11) 

This Wn algebra represents only the chiral part of the algebra of generators of the symmetry, 
which governs the theory. Total algebra is a tensor product of the both holomorphic and 
antiholomorphic algebras Wn Wn. 

Basic objects of conformal Toda field theory are the exponential fields parameterized 
by a (n — 1) component vector parameter a 

14 = e("'^), (1.12) 



^Bellow we consider mainly the case of sphere, in order to avoid the problem with moduli. It is useful 
to choose the metric gab = ^ab everywhere except the north pole {z = oo), where the curvature is located. 
Such a choice prescribes the asymptotic (p = — Qlog |2:| -|- . . . at 2: — > oo. 

^More strictly, it becomes to be invariant under the combined Weyl transformation: gab — » ^{x)gab and 
(fi ^ if — Q\ogQ.{x). 
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w^muzrz) = ^ r T ' + • • • , (1-13) 



which are the spinless primary fields. They have the simple operator product expansion 
(OPE) with the currents W''(4). Namely, 

here . . . means the contribution of less singular terms. Similar OPE's with antiholomorphic 
currents W (^) are also valid. The quantum numbers w^''^ (a) possess the symmetry under 
the action of the Weyl group W of the Lie algebra s[(n) (which is generated by reflections 
in the hyperplanes perpendicular to the simple roots e^) [23] 

w^''\a) = wf\a) = w^''\Q + s{a-Q)), seW. (1.14) 

In particular, 

w^^Ha) = A{a) = ^-^2^Z^ (1.15) 

is the conformal dimension of the field Va- Equation (1.14) suggests the idea, that the 
fields related via the action of the Weyl group should coincide up to a multiplicative factor. 
One of the important properties of TFT is that it is really true 

VQ+Sia-Q) = Rs{a)Va , (1.16) 

where i?s(a) is the refiection amplitude, which was found in [24] 

Rs{a) = A{Q + s(a - Q))/A{a), 

A{a) = (7r/x7(6'))^ nr(l-%-Q>e))r(-6-i(a-Q,e)). ^^-^'^^ 

e>0 

In Eq (1.17) the product goes over all positive roots. 

Multipoint correlation functions of the exponential fields 

{VaAzi,Zi)...Va^izi,Zl)) = J [Vip]e-^^^^ Va^Zl, Zi) . . . Va^{zi, Zl) (1.18) 

are the main objects of the theory. One of the most important problems in TFT is to 
find these quantities. This problem is nontrivial due to the exponential interaction term in 
the Lagrangian (1.1). One can try naively to explore perturbation theory in cosmological 
constant fi. However, pertubatively, correlation functions (1.18) are equal to zero unless 
the "on-shell" condition 

I n-l 

Y,Oij + hY,Skek = '^Q (1-19) 

3=1 k=l 

with some non-negative integers Sk is satisfied. Alternatively, one can perform zero mode 
integration [25]. Namely, let us define a zero mode cpo of the field ip: <p = cpo + ip with 
the condition that J x (f = 0. The integral in Eq (1.18) over the zero mode cpo can be 
transformed to the Euler integral. As a result, after integration we obtain 



{VaAzi,Zl)...Vai{zi,Zl)) 



VaAzi,Zl)...Va,{zi,Zi), (1.20) 
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with 

Sk = ^ . 

Here vectors ujk are the fundamental weights of the Lie algebra sl{n)^. Integration in Eq 
(1.20) is performed in the theory of a free massless {n — 1) component scalar field with the 
action 

5n = ^ y" ida<pfd''x. 



Equation (1.20) has no meaning if all numbers Sk are general. However in the resonance 
situation, when all numbers Sk are non-negative integers, the gamma functions in the right 
hand side of Eq (1.20) have simple poles in each of the variables {2Q — ^aj,ujk) for 
A; = 1, . . . , n — 1 and we can treat the main residue in these poles (the residue in each of 
these poles) as the corresponding free field integrals. Namely 

res ... res (V^^ (zi, zi) . . . Fq.^ (z;, z;)) = 

{2Q-J2oij,uii)=bsi aj,LLin-i)=bs„-i 

= r(l^ai(^l, ^l) . . . VaM^mQir ■ ■ ■ {Qn-lY-'h, (1-21) 

Si- ■ . . Sji—1- 

here (. • • )o means average over the free massless fields. In Eq (1.21) we have introduced 
the notations for the so called screening charges 

Qfc = J e^^^'^^'^''^^^, fc = l,...,n-l. (1.22) 

Correlation function in the r. h. s. of Eq (1.21) can be calculated using the Wick rules 
in the free field theory together with integration over the position of all screening fields 

Equation (1.21), which was obtained from the classical arguments, modifies in quantum 
case. Namely, if the screening conditions 

I 

{2Q-Y,aj,i^k)=bsk + b-^~Sk (1.23) 

are satisfied for any two sets (si, . . . , s„_i) and (si, . . . , Sn-i) of non-negative integers, then 
the correlation function (1.18) admits a pole in each of the variable {2Q — '^a,ujk) with 
the main residue being expressed in terms of free field correlation function 

res ... res (zi, zi) • • • T4( (^/, ^0) = 

(2(3-21, ai,u!i)=bsi+b-^si y^Q-Yl ai,Wn-i)=6s„_i-|-b-is„_i 

^_^^S1+---Sn-1 ^_^^S1 + ...S„-1 

Sl!...S„_l! Si!...Sn_i! 

X {VaAzi,zi) . . . V^,{zi,zi){Qiy^ . . . (Q„-i)^"-UQir • • • (Qn-i)'"-^)o- (1.24) 
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They are defined as a dual basis to the simple roots {ei,u}j) = Ji- 
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In Eq (1.24) we have introduced the notation for the dual screening charges 

Qk = J e^~^^<'k,<Pk)ff^ (1.25) 
and for the dual cosmological constant 



Operators Qk and Qk have an important property, that they commute with all generators 
of the both holomorphic and antiholomorphic W algebras. In this paper we will consider 
for simplicity the case, when all numbers .s^ = 0. It is reasonable to suppose, that the 
screening condition (1.23) defines up to the Weyl transformation (1.14) all possible poles of 
the correlation function (1.18), as a function of the parameters ctfc. One should emphasize, 
that a simple pole in the correlation function (1.18) appears if at least one screening 
condition (1.23) is satisfied. 

Knowledge of two-point and three-point correlation functions of the primary fields 
Va is the first step for the calculation of the higher multipoint correlation functions of 
the theory. In the sl(2) case (Liouville field theory), this knowledge together with the 
statement, that conformal blocks are completely determined by the conformal symmetry, 
allows us, in principle, to compute any multipoint correlation functions in this theory [2]. 
In the s[(n) TFT case for n > 2 the situation is more complicated and we need more data 
(see for example [19]). 

Two-point correlation function in TFT normalized by the condition 

(F«(^)F2Q-a(0)) = |z|-^^(«). (1.27) 

All other non-zero two-point correlation functions can be obtained from this correlation 
function by the Weyl reflection (1.16). For example 

(K«(^)^a*(0))= 1^^, (1.28) 

here R{a) is the maximal refrection amplitude defined as 

A{2Q - a) 

R{a) = (1-29) 

with A{a) given by Eq (1.17) and conjugated vector parameter a* defined as 

(",efe) = (a*,e„_fe). (1.30) 

Much more complicated object - three-point correlation function has standart coordi- 
nate dependence due to the conformal invariance of the theory 

i\i ( - ^T/ ^ - I - w C(ai,a2,a3) 

(Fai {zi,Zi)Va2 [Z2, ^^2) Vag (2:3, ^s)) 



^^2|2(Ai-hA2-A3)|2^3|2(Ai-hA3-A2)|223|2(^2-hA3-Ai)- 

(1.31) 
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All non-trivial information about the operator algebra of the primary fields Va of the model 
is encoded in the constants C(q:i, 02, as). According to Eq (1.21) if the parameters ai, 02 
and as satisfy the screening condition 

q;i + 0:2 + «3 + bsiei H h 6s„_ie„_i = 2Q, 

function C{ai, 0:25 as) will have a pole in each of the variables {2Q — ai — a2 — as,ujk) and 
we can define the main residue in these poles in terms of Coulomb integral"^ 

res ... res C(ai, q;2, as) = 

{2Q-'^ai,uJi)=bsi {2Q-'^ ai,u>n-i)=bs„-i 

= (-7r/x)^i+-+^"-iJ,,...,„_,(ai,a2,a3) (1.32) 

with 



-fsi...s„-i(ai,a2,Q;3) = J dusiih) . . . (i/is„-i(in-i)x 



n— 1 Sfe n—2 



k=l j=l 1=1 

Ai) 

Asismi'th ^m) for I 7^ m are defined as 



here t^^^ is the coordinate of the j'-th screening field e''^'^'"'^^ and quantities 'Ds^{tk) and 



^SlSm \ 

T^^M = n 14^^ - As,sju,tm) = n n if - 4'^ • (1-34) 

i<i' i=H'=l 

Throughout this paper we use the notation for the measure of integration 

dt^sSk) = -^,f{dhf . (1.35) 

In the case of algebra sl(2) Coulomb integral (1.33) is known also as two-dimensional 
generalization of Selberg integral. It can be calculated explicitly in terms of F-functions 
[26, 27, 28] (see also Refs [29, 30]). Unfortunately, it is not clear how to calculate integral 
(ai, a2, as) for arbitrary parameters a^ in the case of general n > 2, but if one of 
the parameters satisfy the special condition, for example 

as = >i:uJn-i, (1.36) 

then the integral (1.33) can be carried out explicitly in terms of F-functions (see appendix 
A). Namely, the integral /si...s„_i (ai, >f^n-i) is non-zero only if si < S2 < • • • < 
In order to write down an answer we define an auxiliary function 

I n 

4 = n ^(-^^') n ^(^(Q - - ^k) - ib^hm - a2, hj - hu) - ib^), 

1=1 j>k 



''in the integral (1.21) we can set using the projective invariance zi = 0, Z2 = 1 and zs = oo. 
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with 



(1.37) 



r(i-x)- 

Integral (1.33) equals in this case (one has to remember, that parameters ai, 02 and x are 
subject to the condition ai + a2 + xa;„_i = 2Q — bsici — • • • — 6s„_ie„_i) 



-fsi...s„_i(ai,Q;2,>^t^n-l) = 
-1 



It is easy to check, that function 



SlH hSn-l 

n 

j=0 



T3S1 T>S2—S\ r>Sn—l—Sn—2 



(1.38) 



C{ai,a2,xujn-i) = TTHjQ) )h 



(2Q-E«,;,p) 



(T(6))-iT(x) n T((Q-ai,e))T((Q-a2,e)) 

e>0 ^ ^ ^ ^ 

nT(f + (ai-Q,/i,) + (a2-Q,/Ji)) 



(1.39) 



which was proposed in [20], satisfies the condition (1.32) at this special case. Here T(a;) is 
the entire selfdual function (with respect to transformation h ^ ^/b), which was defined in 
[17] by the integral representation 



logT 



{x)= [ 
Jo 



T 



b + b-^ 



-X ] e * - 



sinh^ ( ^ 



X] 2 



sinh ^ sinh ^ 



(1.40) 



This function satisfies functional relations 

T(x + b)= ^{bx)b^-^''^T{x), 
T(x + 1/6) = 7(x/6)62^/^-^T(x). 



(1.41) 



and in fact is completely determined by them for the general real values of the parameter 
6 up to a multiplicative constant, which is fixed by the condition 

6 + 5^1 



This function was firstly introduced by Barnes [31], as a generalization of ordinary Gamma 
function and in the semiclassical limit {b 0) it has an asymptotic 



T{by) b^-y 



as 6 — 0. 



(1.42) 



T(6) m 

One can easily check, that the correlation function (1.39) is consistent with the reflection 
identification of the exponential fields (1.16). Due to the symmetry reason, formula similar 
to (1.39), but with as = kooi is also valid^. Note that the condition (1.36) is crucial at 



^One has to change hk ^ hi = —hn+i-k in (1.39). 
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this point and the general formula for the three-point correlation function is much more 
complicated. 

One of the important sets of fields in TFT form so called completely degenerate fields 
[23] . Completely degenerate fields Vq, in TFT are parameterized by two highest weights Jli 
and ^2 of the finite dimensional representations of the Lie algebra sl{n) and correspond to 
the value of the parameter a (up to Weyl transformation (1.14)) 

a = -bQi - li}2- (1-43) 



These fields posses an important property that in their operator product expansion with 
general primary field Va appear only a finite number of primary fields V^/ with their 
descendant fields 

^-6Qi-fe-iQ2^« = XI ^-6^1 -6- 1^2, a 
s,p 

here by square brackets we denote the contribution of the descendant fields and introduce 
the parameter a'gp as 

a'^p = a - bh^' - b'^h^^ . (1.45) 

In Eq (1.45) are the weights of the representation J7 and C'"5Qj_(,-iq2 a denotes the 
structure constant of the operator algebra. During this paper we will consider for simplicity 

the case 0,2 = 0. 

General structure constant of OPE C^f^^j defined as 

Caf,a2 = Ciaua2,2Q-as)=R{a3)C{aua2,al), (1.46) 

where R^a^) is the maximal reflection amplitude given by Eq (1.29) and conjugated pa- 
rameter 03 is defined by Eq (1.30)^. Strictly speaking, structure constant with completely 
degenerate field defined by Eq (1.46) as 

C-mS = Ci-b^i, a,2Q-a + bhf^) (1.47) 
will be infinite because general weight h^^ of the representation fli has a form 

n-l 

h^^=n,-J2sjej (1.48) 
j=i 

with some non-negative integers Sj and hence the sum of all parameters in the three- 
point correlation function in the r. h. s. of Eq (1.47) satisfies the screening condition 

(1.19). In this case one should treat the structure constant C^^q^^^ as the main residue of 
the corresponding three-point correlation function. This residue is given by the Coulomb 
integral (1.33). Namely 

C-bS = (-7r/x)^i+-+^"-iJ,,..,„_,(-60i,a,2Q-a + 6/i^i). (1.49) 

®We remind, that parameters a* and 2Q — a are connected via Weyl transformation (1.29). 



(1.44) 
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The complexity of these structure constants^ depend drastically on the multiplicities of the 
corresponding weights h^^ . 

To illustrate this fact we give here some basic structure constants, which can be ex- 
pressed in terms of known functions. Let us consider, for example, the case Oi = cok 
corresponding to the fc-th fundamental representation of the Lie algebra sl{n). We denote 

(k) 

as the set of weights hi ' of the fundamental representation tt^ with highest weight u>k 
{hi^^ G Hf.). Then the operator product expansion of the field V^hui^, with arbitrary field 
Va due to Eq (1.44) has a form 



(k) 



To describe the structure constants in this expansion we denote as T^.^ the set of 

(k) 

positive roots e such that e + hg € Hk- Then 



^ ^ ^(_62); n^^(i + ;,2 + ^(^_Q^,))- ^'-^'^ 

This result has been derived from the free field integral (1.49) using the same technique, 
which was used in appendix A to derive Eq (1.33). 

Another interesting situation, when the integral (1.33) can be calculated exactly is 
the structure constants with degenerate field K-6eo (^o = X]fc=i is the maximal root 
corresponding to highest weight of adjoint representation). This operators plays a role of 
intcgrable perturbation of the theory, which moves conformal TFT to the massive affine 
TFT. The operator product expansion of the field F-beo with general primary field Va has 
a form 

F_6eoK, = C7^,,„,, [Va] + J2 CXta i^c-be] , (1.52) 

e 

where the sum goes over all roots of s[(n). The diagonal structure constant C^^e^ „ can be 
represented as 

^ h\h ^i-b'h{^+b'+b{a-Q,hj-h.))^^ ^'-''^ 

where functions J-i{a) can be expressed through the higher hypergeometric functions at 
unity „F„_i(l) as 

fr[ fJi (1 + b{Q - a, hj - hi))k 

where 

{x)k = x{x + l)...{x + k-l). (1.55) 



'^We study these structure constants in more details in forthcoming paper [22]. 
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For the positive roots e = hj — hi with i > j the structure constant C^^^^^^ is given by the 
product of 7-functions 



-beo,a yj^_l)2'^ J 
j 

X 



T-r^ -f{b{Q-a,hk-hj)-b'^) y} -r{b{Q - a,hi - hk) - b"^) 
l}^ 7(1 + b{Q -a,hk- hj)) 7(1 + biQ -a,^- %)) ' ^ ' ^ 



While the structure constants for the negative roots can be expressed through the structure 
constants for the positive roots (1.56) as 

Clt^la = R'\^)R{^ + be)C'^:,-l%,, (1.57) 

where R{oi) is the maximal reflection amplitude (1.29) and a' = a + be. 

An important point should be emphasized here. It follows from Eq (1.53), that the 
structure constant C*";,,.,, „ is expressed in terms of higher liypcrgcomctric functions nFn-i 
at unity (1.53), while the structure constants C°i^^^^ have more simple form and are 
expressed in terms of product of 7-functions (1.56). The difference between these two 
cases is related with the fact that field with zero weight in the adjoint representation with 
highest weight eo appears with multiplicity (n — 1), while the weights corresponding to 
the roots e appear with multiplicity equal to 1. The same is true for the fundamental 
representation with the highest weight Wfc, where all weights hf^ of this representation also 
appear with multiplicity 1 and as a result the structure constants (1.51) are expressed in 
terms of 7-functions. We see, that the fact that some weights have multiplicity more than 
one makes the situation more difficult. In the 5l(2) case (Liouville field theory) it does not 
happen because all weights appear with multiplicity one. 



2. Differential equation 

Three-point correlation function (1.39), which was derived in section 1 by the calculation of 
the Coulomb integrals, can be obtained also from rather different arguments. The idea is to 
explore the associativity condition of the operator algebra and to use the special properties 
of degenerate fields. This approach was proposed in Ref [21] in order to find the structure 
constants in the Liouville field theory (5l(2) TFT). Here we will consider in details the case 
of 5t(3) TFT, as the next step of complexity. 

The chiral part of the algebra of symmetries in this case consists of two currents of the 
spin two and three^ 

00 J CO „ , 

W^(.) = r(.)= ^ and W^(.) = I^(.)= ^ (2.1) 

n=— 00 n=— 00 

®This basis of currents is slightly differs from the basis defined by Miura transformation (1.8). Basis 
(1.8) is more convenient, because commutation relations of the W algebra are bilinear. In Eq (2.1) the 
current W is primary field with respect to Virasoro algebra and differs from the corresponding current in 
Eq (1.8) by adding term proportional to T'. 
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The Laurent componets and Wk form closed W3 algebra with the commutation relations 
[3, 32] 

[Ln, Lm] = (n - m)Ln+m + Y2 ^""^ ~ 'n)Sn-m, {2.2a) 
[Ln, Wm] = {2n - m)Wn+m, (2.2b) 

c 16 
[Wn, Wm] = 3751 ("^ - - 4)n5n -m + 22^5^^" ~ m)A„+„ + 

+ {n - m) (^-^{n + m + 2)(n + m + 3) - ^(n + 2)(m + 2)^ Ln+m, (2.2c) 

here 

00 ^ 

k=—oo 

X2i = (! + /)(!- /) X21+1 = (2 + 0(1 - 0- 

This algebra is not Lie algebra due to the quadratic terms in the r. h. s. of Eq (2.2c), 
however, as was noticed by A. Zamolodchikov [3], the Jacoby identities are satisfied. 

The operator product expansions of the holomorphic currents (2.1) with the primary 
fields Vr, has the form 



(2.3) 



here 

{2Q - a,a) 

A(a) = (2.3a) 

is the conformal dimension and 



/ 48 

w{a) = id 22 + 5c ~ ^' ~ ^' ~ ^' {2.3h) 

is the quantum number associated to the W{z) current. Along this section we omit some- 
times (where it is not important) the z dependence of the fields V^. In Eq (2.3) we introduce 
the notations W^iVa{z) and W-2Va{z) for the W descendant fields. Using Eq (2.3) one 
can obtain Ward identities 

N 



N / 

{W{z)V^{z^)...Vm{zn)) = Y.\ 

k=l \ 



\^ ^k) \^ ^k) \^ -^kj j 

(2.4b) 

Let us explain our notations. For example 

t^!?(^i(2i) . . . Vn{z^)) = (Vi(zi) . . . W.xVk{zk) . . . Vn{zn)). 
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One should emphasize, that contrary to the Virasoro generators operators W^k generally 
speaking do not act on correlation functions as some differential operators. This important 
difference explains the essential complication, which appear in the analysis of the s[(n) 
conformal TFT for n > 2. 

The transformation laws for the currents T{z) and W{z) under the holomorphic sub- 
stitution z ^ f{z) have a form ^ 



T(z) 



(2.5) 



The condition, that infinity is a regular point leads to the following asymptotic condition 
for the currents T{z) and W{z) 



1 



T{z)r^-j at z 



oo 



(2.6a) 



and 



W{z) ^ at z ^ oo. 



(2.6b) 



It follows from the asymptotic (2.6a) of the current T(z), that correlation functions of 
the primary fields satisfy certain set of differential equations, which restrict their possible 
coordinate dependence [2]. In the particular cases of two and three points, correlation 
functions are completely determined by them up to a numerical factor. For the case of 
two-point correlation function these differential equations put the limitation on it. Namely, 
two-point correlation function is non-zero only if the dimensions of two fields are equal. 
One can consider the asymptotic (2.6b) of the current ^^(z) in a similar way. Applying 
(2.6b) to the Ward identity (2.4b), we obtain five algebraic equations, which also restrict 
possible form of the correlation functions. Let us illustrate, how does it work in the case of 
two-point correlation function. These five algebraic equations connect different correlation 
functions, which enter in Ward identity (2.4b). Namely, we obtain a system of equations 






Wi + W2 




1 



1 1 \ 

1 Zi Z2 

v2 ^2 



2zi 2Z2 zf Z2 



3{wiZi + W2Z2) Szf 3z| zf z\ 
V ^{wxz\ -h W2zl) 4zf 44 zf z\ ) 



I (Fl(^l) y2{z2)) \ 

(w^_iyi(^i) 1^2(^2)) 

(yi(zi) W.XV2{Z2)) 
(H^-2Vl(zi) 1^2(^2)) 
\{VX{ZX)W.2V2{Z2)) ) 



= 



(2.7) 



This algebraic system has a non-zero solution if the determinant of the matrix above equals 
to zero for any points zi and Z2- A simple calculation leads to 



det = —{wi + W2){zi2)^ 



(2.8) 



^T{z) does not transform like a tensor, but is shifted by the Schwartz derivative, which is defined as 
{f^z} = f" If - 3/2(/"//')^, while W{z) is really a tensor as follows from Eq (2.2b). In the case of s[(n) 
algebra for n > 3 it is also possible to choose currents W''(2;) is such a way, that they will bo primary with 
respect to stress-energy tensor, i. e. will transform like a tensors under the change of variables. 



-14- 



It means that the correlation function {Vi{zi) ^2(^2)) is zero unless wi = —W2- As a result, 
we obtain the following form of the two-point correlation function 



mzuz,)V,{z,,Z2)) ~ ^^f^^ri;-"^ (2.9) 

\Z12\ ^ 

Omitted multiplicative constant in (2.9) depends only on the particular normalization of 
the fields. 

To extract the information about the fusion rules it is reasonable to study completely 
degenerate representations of the W3 algebra (2.2). Namely, if parameters {A{a),w{a)) 
corresponding to the field Va take one of the four values 

A = — ^-1 = — — 1, 2.10a) 

3 27 36+f ^ ^ 

or in terms of parameter a (modulo Weyl transformation (1.14)) 

a = —bcjk or a = —^lo^ k = 1,2. (2-11) 

Then this field exhibits three null-vectors [32, 33, 34] 

XI = (w., - ^^-1) Va = 0, (2.12a) 

/ 16w o 12w ^ ^ 3w; (A-3) ^ \ ^ 

= [^-' ~ A(A + 1)(5A + 1)^-1 + A(5A + 1)^-^^-^ + 2A (5ATT)^-V " ^ 

(2.12c) 

The next natural step is to investigate, how equations (2.12a), (2.12b) and (2.12c) put 
the limitations on the three-point correlation functions, i. e. we want to define the fusion 
rules. Let us consider three-point correlation function {V {z)Vi{zi)V2{z2)) , where field V{z) 
is degenerate field with parameter (2.11) and fields Vi{zi) and V2{z2) are some arbitrary 
primary fields. In the Ward identity (2.4b) for this case participate seven functions: 

{V{z,z)Vi{zi,Zi)V2{z2,Z2)) (2.13) 

and also six functions, which can be obtained by the application of the operators W-i and 
W-2 to the fields V, Vi and V2. Due to conformal invariance, the coordinate dependence 
of the three-point correlation function is known explicitly 

{V{z)Vi{zi)V2{z2)) - (^ - Zi)(^2-Ai-A)(^ _ ^^)(Ai-A2-A)(^^ _ ^^^(Ai+A^-A)^ (2.14) 
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Applying equations (2.12a) and (2.12b) to Eq (2.14) we can express correlation functions 

{W.iV{z)Viizi)V2{z2)) and (1^-2^(^)^1(^1)^2(22)) as 

{W^,Viz)V,iz,)V2iz2)) = {^J^^ + ^J^^) (^(-)^i(-i)^2(.2)) 

(2.15a) 

and 

/w T/^M/A M/^ ^\ /(A + Ai-A2)(A + Ai-A2 + 1) ^ 
{W-2V{z)Vi{zi)V2{z2)) = . . . I ^ r^^ ^ + 



[A(5A + 1)V i^-zi)^ 



^ 2(A + Ai-A2)(A + A2- Ai) ^ (A + A2 - Ai)(A + A2 - Ai + 1) 



{Z-Zi){z-Z2) {Z-Z2 



\2 



6w{A + 1) f{2Ai + A - A2) ^ (2A2 + A - Ai) (Ai + A2 - A) 



A(5A + 1)V (z-zi)^ {Z-Z2Y {z - zi){z - Z2) , 

X {Viz)Vi{zi)V2{z2)). (2.15b) 

Similar to the case of two-point correlation function we obtain five equations, which follow 
from the asymptotic condition (2.6b). The determinant of the corresponding matrix should 
be zero for any points z, z\ and Z2- It gives the equation 

12u;(Ai - A2)^ - 3u;(A + l)(Ai + A2) + A(5A + \)(wx + ^2) - 4wA(A - 1) = (2.16) 

We should take also into a account Eq (2.12c) and put {x'i{z)'V\{z\)V2{z2)) =0. As a 
result, we obtain the second algebraic equation 

32u'(Ai - A2)^ - 12u;(A + 1)(A? - A^) - w{\h^ - 18A - l)(Ai - A2) + 

+ A(A + l)(5A + l)(w;i -u-a) = (2.17) 

The equations (2.16) and (2.17) define the fusion rules in our model (one should fix pa- 
rameters (Ai,tt;i) and find admissible parameters (A2,tt;2)) after that. If we parameterize 
(Ai,i(;i) = (A(a), i(;(q;)) and (A, it;) = (A(— bwi), iti(— ftwi)), then three solutions to the 
equations (2.16) and (2.17) are 

A2 = A(ai - hhj) W2 = -w{ai - bhj) j = 1, 2, 3, (2.18) 

where A(a) and w{a) are defined by Eqs (2.3a) and (2.3b). There are analogous formulae 
for the other completely degenerate fields. We see, that these fusion rules coincide with 
those obtained in [23, 32, 34] and coincide with the fusion rules (1.50) for the Lie algebra 
sl(3) (see section 1). 

Having such rather simple fusion rules (2.18) one can hope that the four-point correla- 
tion function, which contains completely degenerate field, will satisfy differential equation 
of the third order. Unfortunately, this is not the case. Consider, for example, the correla- 
tion function 

{V{Z, z)Va, {Zl, Zi)Va, {z2, Z2)Va, {Z3, Z3)) ■ (2.19) 

Here V{z) is the degenerate field with the parameter a = —huoi. Firstly, one should notice 
that the number of equations in this case is not enough to write down the differential 
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equation. Really: in this case the number of correlation functions in the Ward identity 
(2.4b) is nine. These nine correlation functions satisfy five projective Ward equations plus 
three equations (2.12a), (2.12b) and (2.12c), which arise in the case, when one of the four 
fields is completely degenerate. Total number of equations is eight. Hence, these equations 
allow us only to express all correlation functions in terms of only one correlation function, 
but not to write down the differential equation for this function. Therefore we need at 
least one more additional condition, which connects different correlation functions in Eq 
(2.4b) together. 

Let us suppose that one of the fields V^^, or Va^ in correlation function (2.19) is 
partially degenerate. For example we suppose, that quantum numbers A3 and ws of the 
field satisfy the relation 



9,4 



which can be written, as a condition on the vector parameter as (modulo Weyl transfor- 
mation) 

03 = xwa (2.21) 

with arbitrary coefficient x. Corresponding field V^^a^i^ satisfies the null vector condition at 
the first level 

W-i - ^^-1) V^, = 0. (2.22) 

Under this assumption, correlation function (2.19) satisfies differential equation of the third 
order. In order to write it explicitly, we define function G{x,x) as 

{V{z,z)VaAzi,zi)Va,{z2,Z2)V^,{zs,zs)) ~ |x|2''("i''^i)|l-a;|^-^^^, (2.23) 

I z Z2 1 

with X being the projective invariant of four points x = {l-zl) ^^^^ ~ means that 
we have omitted factors independent on the coordinate z. We derive from Eqs (2.12) and 
(2.22) that function G{x,x) satisfies generalized Pochgamer hypergeometric differential 
equation of the type (3, 2)^° 

X ( x4- + ^1 ] + ^2 ) + A3 

ax \ ax \ ax 



- [x-^ + Bi-l \ \x— + B2-l]x— 



with 



and 



dx ) \ dx J dx 

hx 2 



G{x,x) = fd (2.24) 



Ak='^-^-^h^ + b{ai - Q, hi) + b{a2 - Q, hk), (2.25) 



Bi = l + 6(ai -Q,ei), 

B2 = 1 + b{ai — Q,ei + 62). 



(2.26) 



°Of course, the same differential equation with x being replaced with x is also valid. 
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Three linearly independent solutions to Eq (2.24) with the diagonal monodromy around 
the point x = have a form 



and 



here 



TP I A-1A2 A3 



-B1-B2 



5i(Bi + 1)^2(52 + 1) 



2! 



(2.27a) 



G2(x) = .i-^^F |x) , (2.27b) 

G,{x) = x'-^^F (^'-^^^tk^B^nf-'^''' > (2-27C) 



^ ^ AiA^A, ^ ^ Ai(^i + 1)^2(^2 + 1)^3(^3 + l) x^ , ^2.28) 



is the hypergeometric function of the type (3,2). 

Let us take now into account the antiholomorphic part of the correlation function 
(2.19). We wish our correlation function be invariant with respect to moving point x 
around point 0. The invariant combination G{x, x), which defines the four-point correlation 
function (2.19) has a form 



G{X,X) = X]<^-Lrai^("l ~ ^^J' "2, XUJ2)Gj{x)Gj{x) 



(2.29) 



with C'^l^^'^^^ being the structure constants of the operator algebra. Now we should impose 
the condition that this correlation function remains invariant, if we move point x around 
points 00 and 1. Evidently, it is sufficient to provide this invariancc around point 00, 
because contour surrounding points 1 and 00 can be transformed to contour surrounding 
point 0. 

There is another set of solutions to the equation (2.24), which have diagonal mon- 
odromy around the point x = 00. 



TT ( \ _ -Ai p f Ai l+Ai-Bi I+A1-S2 
Ill[X)-X P\ I+A1-A3 



and 



m(x) - x-^^ F ( 1+A3-B2 
n^yx) — X r \ i_,_^3_^j 1+A3-A 



A3 
-A2 



(2.30a) 



(2.30b) 



(2.30c) 



Of course, these two bases (2.27) and (2.30) of the solutions to Eq (2.24) are linearly 
connected. Using Mellin-Barnes representation for the generalized hypergeometric function 
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one can obtain the relation between them. For example 



r(^i)r(^2)r(^3 



Ai A2 A3 
Bi B2 



.^^ r(^i)r(^2-^i)r(^3-^i) j. [a, i+a,^b, i+a.-b^ 



T{Bi - Ai)T{B2 - Ai] 



I+A1-A2 l+Ai-As 



- I + 

X 



1+A2-B1 A2 1+A2-B2 



1+A2-A1 1+A2-A3 



-A3 



r(Bi - A2)T{B2 - A2) 

TiA,)TiA,-A,)TiA2-A^ ^ f i+As-b, 



X 



+ 



r(Si - A3)r{B2 - As) 



y 1+A3-A1 1+A3-A2 



. (2.31) 



Our correlation function has to be also single valued at the point x = 00. Hence it must 
be represented by the diagonal bilinear form 



G{x,x) = Xl<^-Li,4'^(°'i' "2 - bhj, HUJ2)Hj{x)Hj{x). 



(2.32) 



The necessary conditions of the validity of the both s-channel (2.29) and t-channel (2.32) 
decompositions are 



bin- 02. ^^2) _ nLi -^i-^&i{BL A,} -{1 -B, + B2) 



CXua,C{ai-bh2,a2,xL02) 7(^1)7(52) 7(^1-1) ' 

C-bj' a\C'(ai - bhi, a2, HU2) nLi l{^k)l{B2 - Ak) 7(1 - ^2 + B^) 



(2.33) 



CXf,a\^^(«i-«'^3,a2,xa;2) 



l{BMBi) 



l{B2 - 1) 



Of course, functional equations similar to (2.33) with ai being replaced by a2 are also 
valid. 

One can expect, that differential equation similar to (2.24) will take place in the si{n) 
case too^^. The condition (2.21) undergoes natural modification 

Q3 = xwn-i. (2.34) 

Let us consider correlation function 

(y_5c.i(^,^)v;,(o)y„,(oo)y^^„_,(i)) = |x|2^("i''^i)|i -x|^g(x,x). (2.35) 

Function G(.x, x) satisfies generalized Pochgamer hypergeometric differential equation of 
the type (n, n — 1) in each of the variables x and x 



X ( x-^ + Ai 
ax 



d 

X— h An 

ax 



X— + B1-I 

ax 



X— + Bn-l 

ax 



Uxj- 
ax- 



G{x,x) = (2.36) 



^^We do not give here the strict algebraic proof of this fact for s[(n) with n > 3, but the generaUzation 
is very straightforward. 



- 19 - 



with 



and 



Ak = 



hx (n — 1) 



n 



n 



I? + b{ai - Q, hi) + b{a2 - Q, hk), 



Bk = l + h{ai - Q, ei H h Cfe). 



(2.37) 



(2.38) 



The basis of the solutions to differential equation (2.36) with diagonal monodromy around 
the point x = has a form 



(2.39a) 



for A; > 1, 



(2.39b) 



while the dual basis of the solutions with diagonal monodromy around the point x = oo 
can by represented by the functions 



Hkix) = x-^^F(^^_ 



l+Ak-Bi...Ak...l+Ak-Bn-i 
Ak-Ai...l+Ak-Ak-i l+Ak-Ak+i...l+Ak-An 



(2.40) 



where F 



Al... An 



^ ^ xj is the hypergeometric function of the type (n,n — 1). Four-point 

correlation function (2.35) should be single valued function of the variables x and x. It 
means, that it should be represented simultaneously as 



and as 



G{x,x) =^C'^l^^'^i^C{ai - bhj, as, }€Un-i)Gj{x)Gj{x) 



G{X,X) = Xl'^-Li*,a2^("l' "2 - bhj, XUJn-l)Hj{x)Hj{x), 



(2.41) 



(2.42) 



where functions Gj{x) are given by Eqs (2.39) and functions Hj{x) are given by Eqs (2.40). 
Using formula naturally generalizing Eq (2.31) for n > 3, we can connect two bases Gj{x) 
and Hj{x). As a result, we obtain that the condition of the validity of the both f — and s— 
channel decompositions (2.41) and (2.42) for the correlation function (2.35) has a form 



^-CiC'(ai - hhi, 02, HUn-i) _ XTj^i l{Aj)^{B^_i - Aj) H^^t^ 7(1 + Bj - B^-i) 



l{Bk-i - 1) 



(2.43) 



where k = 2, . . . ,n. 

The structure constants C'^^^^'^'^^^ admit the free-field representation [35] 



^ai—bhk 
—boJi,ai 



„ fe-1 

-ai+&fefc(oo) llVbe,iZi,Zi)(fzi)o. (2.44) 
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The expectation value in Eq (2.44) is taken using the Wick rules in the theory of a free 
massless scalar field. This integral, as was pointed out in section 1, can be calculated 
explicitly. The answer follows from Eq (1.51) 

f^ai-bhk _ f TT/X y~^TT ^ {b{ai - Q , kj - kk)) , . 

-ft-i.ai V 7M')y l\j{l + b'' + b{ai-Q,hi-hk))' ^' ' 

Therefore, we obtain from Eqs (2.43) and (2.45) the system of (n — 1) functional relations. 
There is another dual set of functional relations with parameter b being replaced with b~^ 
and cosmological constant n being replaced with dual cosmological constant /x defined by 
Eq (1.26). If parameter b^ real and irrational, then the solution to the both systems of 
equations is unique up to a multiplicative constant, which depends only on the parameter 
X. It is easy to check that proposed in section 1 three-point correlation function (1.39) 
satisfies both of these systems of equations. 

In conclusion of this section, we present the exact expression for the four-point cor- 
relation function (2.35). This correlation function can be expressed in terms of Coulomb 
integral 

(y_6a;i {X, X)V^, {0)Va, {oo)V^^_, (1)) = 



n-1 (2Q 

e>0 



7r/x7(6')6'-'''' 



^ (T(6))--^ T(x) n T((Q-ai,e))T((Q-a2,e)) 



n -f (ai - Q, hi) + {a2 - Q, hj) - bSij) 

ij 

/n—l 
n d^tk \tk\^^'''-'''% - - (2.46) 



k=l 

where tn = ^, o: = —buji + ai + 02 + xcj^-i and parameters and are given by Eqs 
(2.37)-(2.38). This expression for the correlation function can be derived by the analytical 
continuation to the non-integer values of numbers in Eq (1.21)^^, which permits also to 
find expressions for more general correlation functions (see also [29, 30]). In principle, it is 
possible to write down explicit expressions for the correlation functions 

{V-mbo., {X, x)Va, {0)Va, {oo)V^^„_, (1)) 

in terms of finite dimensional integral for m > 1, but the result will have more tedious 
form. We plan to consider these and more general correlation functions in the forthcoming 
paper [22]. 

If we consider the operator product expansion of the field V-hw^ with the field V^i^^_-^ 
in the correlation function (2.46), we find that the coefficient before singularity (1 — x^^l"^ 
defines the three-point correlation function C(ai, 0:2, x.u)n-\ — bui), which is given by the 



^■^It can be also proved, that integral in Eq (2.46) satisfies holomorphic (and antiholomorphic) differential 
equation (2.36). 
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expression 



C{ai, a2, - bui) = 



ra-l 



(2Q-c,p) 



(T(6)r-iT(x) n T((g-ai,e))T((g-a2,e)) 

e>0 ^ ^ 

n T + (ai - g, + {a2 - Q, hj) + bSij) 

ij 

J n - - ir'^^ (2.47) 



where = 1- Integral in Eq (2.47) can be calculated in many different ways. The simplest 
one is to combine Eqs (2.35), (2.46) and (2.41) and express it in terms of hypergeometric 
functions of the type (n,n — 1). As a result we obtain that 



/n— 1 
k=l 



7(^1 



n-l 







-Bfc+Ai ... 1-Bk+A. 
+Bi ...,2— Bfc,... 1— Bfc+B, 



fe=l 



n-l ) 



where 



^/ ^1 \ def 7(-4l)..-7(^nj p/ Ai ...An l-iV 

"^l^Bi ...B„_i; 7(Bi)...7(5„_i) V^i • Bn-ilV 



(2.48) 



(2.49) 



As we see from the above, the four-point correlation function in sl{n) TFT, which 
contains completely degenerate field, satisfies differential equation, only if at least one of 
the other fields is special. Namely, if field Va^ is degenerate at the first level (parameter 
as takes the special value = xa;„_i), then the four-point correlation function satisfies 
FTichsian differential equation of the order n, which can be reduced to the generalized 
Pochgamcr differential equation (2.36) and, hence, it can be represented by the Coulomb 
integral (2.46). Without such a condition, the four-point correlation function seems to be 
more complicated object. One can prove, that for n > 2 it can not be a solution to the 
ordinary Fuchsian differential equation of the order n [20]. However, if the field Va^ is 
degenerate (but not completely degenerate) at the higher levels A; = l,...,n — 2, 

i. e. parameter as takes the special values 



n-2 



as = xun-i - X 'rn-kb^^k 



(2.50) 



k=l 



with non-negative integers m^, then the four-point correlation function can be represented 
by the Coulomb integral of the finite order 



^^We will show it for the case of 51(3) TFT in Ref [22]. 
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3. Classical limit (heavy exponential fields) 



In this section we consider the semi-classical limit b of the conformal TFT. Let us 
define classical field as 

(p = hip. (3.1) 
Its dynamics is described by the classical action 

ra-l 

(9(/))2 + 87r/x62^e('='='*) . (3.2) 
fc=i 



^class 



87r62 / 



In this limit the leading asymptotic of the correlation functions (saddle point asymptotic) 
is governed by the classical action calculated on some specific solution to the equations of 
motion, which follow from the action (3.2). 

We will consider here the case of sl{3) TFT, as an example {sl{2) case corresponding 
to Liouville field theory was considered in [17]), which is already non-trivial. The main 
asymptotic at 6 —> of the correlation functions {Va^{zi,zi)...Vaj^{zN,ZN)) of heavy 
operators with parameters 

ak = 'f (3.3) 

is given by the regularized action^^ 

{Vai {zi,Zi) . . . Vaj^{zN, Zn)) ~ ^'^vi^-^dass ['^('^1 " " " W^l, ^1 • ■ ■ ^AT, Zn)]^ , (3.4) 

here ^(tji . . . r/jvl-^ij zi. . . zn, zn) - real single-valued solution to the Toda equation 

dd(l) = TTfib'^ (^eie^^i-'^) + 626^"^^'^'^^ (3.5a) 

with the asymptotic conditions 

(f) = — 4plog |z| -I- . . . at \z\ — 00, (3.5b) 
(f) = —2r]j log \z — Zj \ + Xj + . . . at z ^ zj. (3.5c) 

In Eq (3.5c) Xj is a. z independent term. The solution to the boundary problem (3.5) with 
positive cosmological constant /j, exists if 

N 

J](%,a;fe)-2>0 A; = 1,2. (3.6) 

i=l 

The regularized action S^^^^^ on this solution can be calculated as follows [17]. By definition 
of the classical regularized action its differential is related with parameters Xj defined by 
Eq (3.5c) in a simple way 

N 



^^The divergences arise from the vicinity of sources corresponding to the insertion of the operators Va 
To obtain finite answer one should re-normahze them. See [17, 36] for the regularization prescription. 
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The constant of integration in Eq (3.7) can be fixed by the condition 



15 



nreg 
class 



N 



In the case of sl{3) TFT it is convenient to introduce the projection of the field (p on the 
fundamental weights Uk, k = 1,2: 

$fe = (<^,a;fc). (3.9) 
In terms of fields equation (3.5a) has a form 

9a$2 = 7r/x62e2*2-*i. (3.10b) 

General solution to the system of equations (3.10) can be obtained by introducing the 
holomorphic currents 

T = {d^i f + {8^2? - d^id^2 - d'^^i - a2$2 (3.11) 

and 

W = (d<^i{d^2f + <9$i<9^$i - ^d^id'^^2 - ^d^^i 

d^2id^if + d^2d'^^2 - ^d^2d'^^i - ^d^^2^ ■ (3-12) 

Using Eq (3.10), one can easily verify that dT = dVM = 0. In a similar way, if we change 
d ^ d in (3.11) and (3.12), we obtain anti-holomorphic currents T and W. It follows from 
the explicit form of the currents T and W, that field e~*i satisfies both holomorphic and 
anti-holomorphic linear differential equations of the third order 

(^-d^ + ^dJ + Jd + \N^e-^' =0, (3.13a) 

-a^ + ^OT + ta + W) e-*i =0. (3.13b) 



Similar equations for e"*^ with changed sign before W and W are also valid-*^^. Differential 
equations (3.13) will play an important role in the following. 

Prom the other hand, equations (3.13a) and (3.13b), being viewed as a system of 
linear holomorphic and anti-holomorphic differential equations with arbitrary functions 
T{z), T{z), W(z) and W(z) can be used to solve the system (3.10). Namely, let = ^kiz) 
are three linearly independent solutions to Eq (3.13a) and = ^k{^) ^^re three linearly 



^''In quantum case this condition means, that correlation function {Va^ {zi, zi) . . . Vq,„ {zn, zn) is trivial 
in the case, then '^ctk = 2(5. Namely, correlation function has a multiple pole under this condition with 
residue expressed in terms of free field correlation function without screening fields (see Eq (1.21)). 

^®It is evident because current T{z) is symmetric and current \N{z) is antisymmetric under the substitution 
1^2. 
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independent solutions to Eq (3.13b)^^. Then we can express the field e as a bilinear 
combination 



3 



Y^^k^k. (3.14) 



k=l 

After that we find the field e"*^ from the equation (3.10a) 

3 

e-*2 = _(vr^^,2^-i (^i9*j- - ^jd'^i){^id^j - ^jd^i). (3.15) 
i<j=i 

Second equation (3.10b) is satisfied only if 

W [^1, *2, *3] W [^1, ^2, ^3] = - (tt/x^^)^ , (3.16) 

here W[^'i,*2,*3] is Wronskian. Henceworth, the solution to the system (3.10) can be 
build up from the solutions of any pair of holomorphic and anti-holomorphic linear differ- 
ential equations of the third order with the condition (3.16). 

In our case, we should solve Eq (3.5a) with boundary conditions (3.5b) and (3.5c). It 
puts the limitations on the possible form of the currents T(z) and \N{z). As a consequence 
of Eq (3.5b), the currents J{z) and W(2;) have asymptotic at infinity 

J{z) ~ ^ W(z) ~ ^ at z^oo (3.17) 
and due to (3.5c) are in fact the rational functions 

N 



k=l 
N 



{Z - ZkY {Z - Zk] 



(3.18) 



.«,/ ^ sr^ I "^k . Dk Ek 



k=l ^ ~ ^'^^^ ~ ^'^^^ ~ 

here parameters 5k and are expressed in terms of vector parameters rjk as 

k 2 ^'^'^> (3.19) 
Wjk = ((%, ui) - {rjk, UJ2)) {{Vk, i^i) - 1) ((%, (^2) - 1) 

and coincide up to a sign with semiclassical limit of the quantum numbers (2.3a) and 
(2.3b). The parameters Ck, Dk and Ek are not defined from the main asymptotic (3.5c) at 
z ^ Zk, but contain information about next subleading terms. In fact, they are not linearly 
independent, but satisfy linear algebraic relations, which follow from the asymptotic (3.17) 
(analog of Ward identities (2.4a) and (2.4b) in quantum case). 

First interesting case is the case of three singular points, which corresponds to the 
semiclassical limit of the three-point correlation function. Let us consider it in more details. 
In this case the number of equations, which follow from the asymptotic of the current T{z), 



Generally speaking functions V&fe and "ifk do not complex conjugated to each other. 
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is enough to find parameters C^- Really, we have three parameters Ci, C2 and and three 
conditions, which follow from the asymptotic of the current T(2;) at infinity. Unfortunately, 
this is not true for the asymptotic of the current W. In this case, we have six parameters -Dfe 
and and only five equations, which appear from the asymptotic W(z) ~ ^. Therefore 
one parameter remains free. Evidently, it corresponds to the possibility to add to the 
current \N{z) the term 



{Z - Zlf{z - Z2Y{Z - ZzY 

with arbitrary coefficient. In order to emphasize this one-parameter freedom, let us fix first 
non-vanishing term of the asymptotic of the current W(2;) at infinity as 

^{z) = TTfi Wi2;i2Zl3(-Zl2 + 2:13) + W2Z21-223(-2;21 + Z23) + 



2^6 L 



+ W3Z3i2;32(2;3i + Z32) + 2AZ12Z13Z23 + O j . (3.20) 



The parameter A, which we call accessory parameter, is not known a priori. Here we arrive 

at the main difFcrcncc with s[(2) case, where the accessory parameters do not appear in 
the case of three singular points [37]. This difference explains at the classical level why the 
three-point correlation function is much more complicated object in higher Toda systems. 
As we will show below, the parameter A can be found, in principle, from rather different 
arguments, which resemble the conformal bootstrap program. 

Using projective invariance of Eqs (3.13)^^ one can rewrite them through the invariants 
of four points 

{Z - Zl){z2 - Z3) _ {Z- Zl){z2 - Z3) 

~ — and 



{z - Z3){Z2 - Zi) {z - Z3){Z2 - Zi) 

For example, Eq (3.13a) will have a form 



(-^' + ^d^Jix) + Jix)d, + W(x)) *(x) = (3.21) 



with 

T(x) - — + _^ S3 - Si- 82 

x"^ {x — 1)2 x{x — 1) ' 



W(X) = ^ + + ^(W, + W2 +W3)fl - ^-^) + - 



x'^ {x — \)^ 2 (z — 1)2/ a;2(a; — 1) 



(3.22) 



Equation (3.21) is the most general Fuchsian differential equation of the third order with 
three singular points 0, 1 and 00 modulo "gauge" transformation ^i{x) —>■ x"(x — l)^"^{x). 



^One can show, that differential equation 



is invariant under the substitution z w{z), ^(z) (^) T{z) -> (^)^T(w) - 2{w,z} and 

~* (^)^^(^)' where {w,z} is Schwartz derivative. 
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The "gauge" is fixed by the condition, that term with second derivative d^'^{x) is absent 
in Eq (3.21). 

In order to solve the problem (3.5), one should find real single valued solution to Eqs 
(3.13). The last requirement is not trivial, because the general solution to Eqs (3.13) does 
not satisfy this property. Let ipk be the basis of the solutions to Eq (3.13a) with diagonal 
monodromy around point z = zi 

^JJ, = {z-z,f+^'^'-P'^''Hl + 0{z-Zl)) k = 1,2,3 (3.23) 

If we write a diagonal bilinear combination 

e-*i = AilV-iP + A2|V'2|' + AalV-al', (3.24) 

such a solution is evidently invariant if we move point z around point zi (contour Ci on 
figure 1). But we need also such an invariance around points Z2 and z^ (contour C2 and C3 
on figure 1 respectively). Let Xk be the basis of the solutions to Eq (3.13a) with diagonal 
monodromy around point z = Z2 

Xk = {z-Z2)^+^'''-P''''^Hl + 0{z-Z2)) k = 1,2,3. (3.25) 

The following formula also should be valid 

e-*i=Ai|xiP + A2|x2p + A3|x3p (3.26) 

with some other constants A^. If the solution e"*^ can be represented simultaneously 
as (3.24) and as (3.26) it becomes single-valued on a total sphere, because the contour 
surrounding point z^ can be transformed to the contour surrounding points zi and Z2, as 
guaranteed by the condition (3.17). As functions V'fc and Xk satisfy the same differential 
equation, they are linearly connected 

V'i = Mijxj. (3.27) 
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Entries of the matrix My are believed to be meromorphic functions of the parameters 6^, 
Wk and the accessory parameter A (for the real values of the parameters 6k, and A 
matrix Mij is real). If we substitute relation (3.27) into Eq (3.24), we obtain unwanted 
cross terms like X1X2 destroying the property (3.26), which guarantees that solution is 
single- valued. So, one should set all coefficients before such terms equal to zero. As a 
result, we arrive to the system of equations 

/M11M12 M11M13 M12M13' 
M21M22 M21M23 M22M23 I I A2 I = 0. (3.28) 

yM3iM32 M31M33 M32M33^ 

The determinant of the corresponding matrix should be zero 




det 



/M11M12 M11M13 Mi2Mis\ 

M21M22 M21M23 M22M23 

\^M3iM32 M31M33 M32M33y 



= 0. (3.29) 



The condition (3.29) can be viewed as an equation on accessory parameter A. Each acces- 
sory parameter A, which solves equation (3.29) gives single-valued solution to the boundary 

Toda problem (3.5). 

Let us try to find a solution to this equation in a special situation, which corresponds 
to the classical limit of three-point correlation function (1.39). Namely, we suppose that 

773 = KUJ2- (3.30) 

In this case one can guess accessory parameter A, which solves equation (3.29): 

A = - {Si -62) -I (wi - W2) (3.31) 

from the rather simple reasoning. The logic is the following: if the condition (3.30) is 
satisfied, the equations (3.13) have the same behavior near the singular points as a hyper- 
geometric equation of the type (3, 2), but do not coincide with it if parameter A is general. 
So, we select such special value of the parameter A defined by Eq (3.31), that these equa- 
tions are identical. One can easily check that the equation (3.29) is satisfied in this case. 
The solution to the boundary problem (3.5) can be obtained in this case in a simple way: 
expression for the field e~*i can be derived by semiclassical limit of four-point correlation 
function (2.46) for n = 3, while expression for the field e~*2 can be obtained from Eq 
(3.10a). Both of them can be expressed in terms of Coulomb integrals (for simplicity we 
set zi = 0, Z2 = 00, Z3 = 1 and z = x) 



2 = £|x|2(''i''^2)|^_l|2-f j ^2^^2y |i_^|2bi|^_y|2b2|y_;^|2b3 |t|2ai|y|2a2^ 



(3.32) 
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with 

^fe = -| + ivi - p, K+i) + {m - p, K), bfc = -1 + ^ - - K) - im - p, K), 



= -1 + I + - P, hk+i) + iV2 - P, hk), bjfc = -| - (?7i - p, hk) - {ri2 - p, h), 

(3.33) 



where h'^ = — ^4_fc and 



^ ^ Pb'^ nLi 7(§ + iVi - P, hk) + {V2 - P, hk)) 

nil UU bil + - hi) + - P, h^))] ' ' 



(3.34) 



^ _ ph^ nil n ■=! [7(§ + (t?i - P, /^i) + (r?2 - P, Z^,))] ^ 

^ nLi7(f + (m-p,%) + (??2-p,%)) 

This solution can be also written through the hypergeometric function of the type (3,2). 
The regularized classical action on this solution -S^^^g can be easily found using Eqs (3.7) 
and (3.8) and has a form 



SlL = iiVi +V2,p) + K- 4) log(7r/x62) + F{k) + ^ F{{p - e)) + ^ F((p - r,2, e))- 

e>0 e>0 

-T.^{^ + im-p, hi) + (r/2 - P, hi)) - F\0) (3.35) 

ij 

with ^ 

F(a;) = /" log7(i)dt. (3.36) 

We see that the classical limit of the proposed three-point correlation function (1.39) is in 

complete agreement with expression (3.35). 

We have found a sohition to the boundary Toda problem (3.5) and an explicit expres- 
sion for the accessory parameter A in the case of five-parametric family (vector param- 
eter rjs restricted by the condition 773 = kol!2)- In general case boundary problem (3.5) 
is rather complicated, because it corresponds to the most general differential equation 
of the third order with three singular points of the Fuchsian type (by projective invari- 
ance this equation can be transformed to equation (3.21)). It is interesting to notice, 
that to the same type belongs differential equation for the four-point correlation func- 
tion {V-b{z)Va^{zi)Va2{z2)Va3iz3)) in the Liouville field theory [2] (sl(2) TFT)!^. It can 
be transformed to the differential equation (3.21) with parameters (one should take into 
account projective invar iance) 

dk = -'^b\6AL{ak)-S-2b^), \^k = -^b^{3(ib^^L{ak)-8b^-m^-9), A = 0. (3.37) 



'^^Here we use standart for the Liouville field theory (LFT) notations, which differ from those used in this 
paper. Namely central charge in LFT equals cl = 1 + 6(6 + b~^)'^ and exponential fields Vc, have conformal 
dimensions AlIq) = a{b + — ot). One should emphasize, that parameter 6 here is formal parameter, 
which does not goes to zero. 
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If we introduce auxiliary parameter g = —b^, then the numbers and vj^ are subject the 
condition 

4 - ^ ^ ^ - 1. (3.3S) 

which can be parameterized in terms of vector parameter r)k, which enter in Eq (3.19), as 

rjk = XkLOi + (Afc - 2g)u;2 k = l,2, 3, (3.39) 

where uji and lj2 are the fundamental weights of the Lie algebra s[(3) and are auxiliary 
scalar parameters. Simultaneous single valued solution to equation (3.21) and to corre- 
sponding antiholomorphic equation can be written in terms of Coulomb integral [29, 30] 
(this solution coinsides up to multiplicative constant with the field e~*i) 



S)|arp^i-^^/3|x-lp^^-4^/=^ f l[\tkf^\tk-lf^\tk-x\^''vl'{t)dHidH2, (3.40) 

k=l 

where 1*2 (i) = |ii — i2p and 



A = -Xi-C, B = -X2-C, C = 1 + 5- ^(Ai + As + Aa). (3.41) 

The dual solution to the differential equation (3.21) with changed sign before current W{x) 
(this solution corresponds to the field e"^^) has a form 



g-<I>2 



2)/|^|2Ai-8p/3|^_l|2A.-8,/3 f Yl\tkf^'\tk-lf'''\tk-xf'''vl''{t)d\dH2 (3.42) 



with 

A' = A + g, B' = B + g, C = C + g, g' = -g. 

Functions e^*i and e^*^ given by Eqs (3.40) and (3.42) define modulo numerical factors 
the solution to the Toda boundary problem (3.5) with three singular points 0, 1 and oo and 
parameters ijk given by Eq (3.39). We do not give here explicit expressions for the numerical 
constants D and D' before integrals (3.40) and (3.42) and for the action calculated on this 
solution, because in this paper wc do not suppose to quantize it. 

The results of this section show, that the boundary problem (3.5) for the s[(3) Toda 
equation is much more complicated than the corresponding boundary problem for Liou- 
ville equation. Even in the case of three singular points one should deal with accessory 
parameters. Boundary problem (3.5) can be reduced to the problem to finding single val- 
ued solution to the holomorphic and antiholomorphic Fuchsian differential equations of the 
third order with given behavior near the singular points Wc have shown that it can be 
reduced to the problem of finding values of accessory parameter A which solves equation 
(3.29). We suppose, that in the domain (3.6) the solution to this equation is unique. 

An interesting question how to find parameter A. We have found it in two different 
cases. First case is when one of the parameters rjk is proportional to the fundamental 



°The most general such equation with three singular points 0, 1 and oo is given by Eq (3.21). 
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weight (for example 773 = nuj2)- In this case accessory parameter A is given by Eq (3.31). 
Another interesting case, which corresponds to the differential equation for the quantum 
field V-b in Liouville field theory gives the value of the parameter A = 0. One has to notice, 
that in both cases the solution to Eq (3.21) is given in terms of Coulomb integrals over a 
plane (Eqs (3.32) and (3.40) respectively), so these solutions are evidently single valued. 
An important problem remains unsolved: how to find parameter A in general case? It is 
difficult to expect, that solution to Eq (3.21) in general case can be expressed in terms of 
finite dimensional integral. Because of that we do not have a efficient procedure to find 
matrix Mij defined by Eq (3.27). We suppose to develop the effective numerical method 
to solve this problem in a future publication. 

4. Classical limit (light exponential fields) 

In this section we consider the semiclassical limit of sl{S) TFT in the opposite case of light 
exponential fields V^^ with parameters 

ak = brik- (4.1) 

The solution to the Toda equation (3.5a) with positive cosmological constant fi in this case 
does not exist, because the condition (3.6) does not satisfied. In order to have a solution, 
it is useful to perform analytical continuation fi — —fi'^^. The leading asymptotic behavior 
of correlation functions at 6 is now governed by the solution to the Toda equation 
with the opposite sign in the r. h. s. 

dd(l) = -TTub'^ (^eie^^i'-^) + e2e(^2,0)^ _ (4 2) 

It is evident, that light exponential fields V^^^, do not affect on dynamics, it means, that in 
this case one has to set 

T = W = T = \A/ = (4.3) 

in Eqs (3.13). General solution to Eq (4.2) expressed in terms of solutions to Eqs (3.13), 
which in this case are the polynomials of degree 2. It is convenient to parameterize these 
polynomials by nine complex parameters ak, hk and Cfe (A; = 1, 2, 3) as follows 

Pk = ak + hz + y 2:^. (4.4) 

The solution to the Toda equation (4.2) is given by 

4P{z, z) = -plogiTTfib") + ei$?(z, z) + e2$^(z, z) (4.5) 

with p being the Weyl vector and 

$0(Z,Z) = - log + |p2p + \P3\^) , ^^2{^,Z) = - log (Ifiip + |^i2|2 + l^g^) , (4.6) 

where the dual polynomials pi are defined as 

Pi = P2P3 - P3P2, P2 = P1P3 - P3P1, P3 = P1P2 - P2Pl- (4.7) 

^'^Alternatively, one can consider the correlation functions with the fixed "area". See Ref [17] for details. 
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Due to Eq (3.16) nine complex parameters a^, and are subject to the SL{3,C) 
constraint 

^ai bi ci^ 

det 02 62 C2 = 1. (4.8) 
\a3 h C3J 

The difference with the semiclassical hmit of the correlation function of "heavy" exponen- 
tial fields considered in section 3 is that the saddle point now is not unique. The action 
is minimazed on any function (p^{z,z) defined by Eq (4.5). So, in order to obtain the 
semiclassical expression for the correlation function of the light operators Vb^^. one should 
integrate over the space of all polynomials (4.4) restricted by the condition (4.8). Namely, 
the semiclassical limit of the AT-point correlation function is given by the integral 

1 f ^ 

l^{Vbr,dzi,Zi)...Vbr,A^N,ZN)) / TT ^^'"'^"^'"''^^^^K, Cfe), (4.9) 

Zo 6-0 J 

where Zq is TFT partition function and dQ{ak,bk,Ck) is the SL{3,C) invariant measure. 

It is evident from Eq (4.6), that fields ^i{z,z), <I>2(-2, 2) and the integral (4.9) have 
SU{3) invariance. Hence, the integral (4.9) is just proportional to the volume of SU{3). 
We use Iwasawa decomposition for the SL{3, C) to fix the gauge. Namely, each SL{3, C) 
matrix can be represented, as a product of SU{3) matrix and uppertriangle matrix with 
unit determinant: 

^ Q a b^ 

SL{3,C) = SU{3) X Que (4.10) 

^0 Tj 

here a, b and c are the complex numbers and g, v and r are the real numbers with the 
condition qvt = 1. In this gauge polynomials and defined by Eqs (4.4) and (4.7) will 
have the following form 



Pi 



gz 



+ az + b, 



P-i 

P3 



gvz 



Pi = VT, p2 = T{gz + a), 

The measure will transform to 

dO,{ak, bk, Ck) = dQ d^a d^b d^c g^dg udu. 



(4.11) 



+ cgz + ac — bv. 



(4.12) 



where dQ is SU{3) invariant measure. 

In the case of three-point correlation function it is convenient to introduce the notations 



(771, Cj) = Xj-, {rj2,ej) = Kj-, (r/3, ej) = aj; j = 1, 2. 
The semiclassical limit of three-point correlation function has a form 



(4.13) 



-^{Vbr,Azi,Zi)Vbr,2{z2,Z2)Vbr,:,{z3,Z3)) (TT/xfe^) 



-Al— A2 — Kl— K2— Cl— <T2 



X J{Xl, X2; Kl, K2; CTi, a2\zi, Zl, Z2, Z2, Z3, Z3), (4.14) 
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with 



J{Xl,X2;i^l,K2;cri,a2\zi,Zi,Z2,Z2,Z3,Z3) = 

= / ^""P {Y1 ^i) + i^k^ki^2,Z2) + cTk^ki^s, 23)) j d^a (fb dPc g^dg udu, 

(4.15) 

where functions ^\{z) are given by Eq (4.6) with polynomials and pk defined by Eqs 
(4.11). The coordinate dependence of the integral (4.15) is fixed by the projective invariance 
and we can set = 0, 2:2 = 1 and z^ = 00. More exactly we consider the limit of the 
integral (4.15) 

5(Ai, A2;ki,K2;o-i,o-2) = lim \z2,\^''''^^''^^ J{\i,\2;i^i,i^2;cri,(^2\^A,zz). (4.16) 

Z3— »0O 

Function 5(Ai, A2; ki, K2; ci, (T2), which defines the semiclassical limit of three-point corre- 
lation function, will be the main object of this section. For convenience, it is better to 
renormalize parameters a, 6, c, g and u 

a — ^ va; b — > rb: c — ^ rc, 

(4.17) 

g — ^ Tg; u ^ tv 

and take into account the condition rvg = 1. Then function CJ(Ai, A2; ki, K2; (^1-1(^2) can be 
represented by the eight-dimensional integral 



a(Ai,A2;Ki,K2;(Ti,(T2) = 4^i+'^2 / — 



dgdv 2 ,2.-2 
d ad bd c 



with 



Zi = l + \bf + \cf, 



1 -I- |a|^ + \ac - 6|^, 



i + \c+iyr + 



b + i'a + 



^4 = 1 + 



a + 



+ 



- H \-ac-b 

2 V 



(4.18) 



(4.18a) 



and 



<^ = o (^1 + 2-^2 + Kl + 2K2 - 2ai - £72) , 
o 

A = ^ (Al - A2 Kl - AC2 + 0"l - (72) . 



(4.18b) 



After non-trivial transformations (see appendix B for details), integral (4.18) can be re- 
duced to the three dimensional Barnes like integral 

a(Al, A2; Kl, «2; CTi, (72) = 4^l+«l+^l-^X 

r(Ai + Ki + (71 - A - 2)r(A2 + /C2 + (72 + A - 2) 

r(Ai)r(A2)r(Ai + A2 - i)r(Aci)r(AC2)r(Ki + ^2 - i)r((7i)r((72)r((7i + (72 - 1) 

X — ^ / dTX dy 4"" r(y)r(s)r(n)r((7i - ■u - s)r(A2 + A - u - s)r((7i - AC2 - A + y)- 

r(c7i -K2 + A1-A-U - s)r(«;i + K2 - 1 - y)r(Ki + Ai - a - 1 - y)r(K2 + (T2 + a - l - y)- 

r(A. - + A + ,) ^'7 -f + ''''f + ^ - 7 + ')^<" - ^) , (4.19, 
r(s + y)r(Ai + Ki+c7i-A-l-n-s-j/) 
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where the integral over variables u, s and y goes along imaginary axis. Integral (4.19) is 
convergent in the domain 

l + im- P,hi) + {r]2- P,hj) + i'n3- P,hk) >0 if {hi + hj + hk, p) > -I, 
^ + im - P, hi) + {112 - p, hj) + {rjs - p, hk) < if {hi + hj + hk, p) < -I, 

where vector parameters r)i, r]2 and r/3 are related with parameters Afc, and ak by Eq 
(4.13). By definition, the integral (4.9) and hence the integral (4.19) should be symmetric 
with respect to substitution X ^ k ^ a and also with respect to substitution 1^2, but 
these symmetries are not evident from it explicit form. We were not able to represent this 
integral in terms of the finite sum of the known functions. Another alternative expression 
for the function 5(Ai, A2; /ci, K2; ci, (72) in terms of three-dimensional integral of Tricomi 
functions is given in the appendix B. Wc sec, that the integral (4.19), which is semi- 
classical limit of the three-point correlation function, is already rather nontrivial object. 
Henceworth, it is difficult to expect that quantum expression will have a simple form. 

Integral in Eq (4.19) simplifies drastically if one of the parameters A^, Kk or ak equals 
to zero. For example, let us consider the limit ai 0. We see, that due to the factor 
r(s)r(n)r((Ti — u — s) in the integrand in (4.19) in this case integral (4.19) develops a pole, 
when we pinch points u and s near the point u = s = 0. This pole cancelles with a zero 
coming from function T{ai)^^ in the prefactor of (4.19) and remaining integral over the 
variable y can be performed using the first Barnes lemma (D.6). As a result we obtain a 
simple expression for the function (4.19) 



5(Al, A2; Ki,K2-A ^2) = A'i^+i^^+^^'P^X 



r(^2) Hoo r(i + {m - p, e))T{l + {^ -p,e))' 
where vectors rik are given by Eq (4.13) and 

e-^. = I ^(f + (''i - + - p, hj)) if {hi + hj,p) > -1 



(4.21) 



r(l - f - (r/i - p, hi) - (r/2 - p, hj)) if {hi + hj,p) ^ -1 



The result (4.21) agrees with the corresponding limit of the three-point correlation function 
(1.39), where ai = br]i, 02 = br]2 and x = ba2- If all vectors rjk are general, the integral 
(4.19) is much more involved object and we plan to study its analytical properties in a 
future publication. 

Another interesting point, where the integral (4.19) can be calculated exactly, is defined 
by the condition ai = —m with integer m, i. e. 773 = iT2a;2 — mui'^'^. In this case function 
(4.19) can be given by triple finite sum, which in general contains 

N{m) = im + l){m + 2){m + 3) 



^^In quantum case field Va with parameter a = a2Ui2 — mbuji is partially degenerate (it has a null-vector 
at the level m + 1). The three-point correlation function with such a field can be expressed in terms of 
finite dimensional Coulomb integral, as it will be shown in Ref [22], and generalizes the answer (2.47) for 
the case of Lie algebra 5l(3). 



-34- 



terms and has a form 



a(Ai , A2; Ki , K2; -m, (72) = 4^i+''i"^r(Ai + Ki - 2 - m - A)r(-m - A)r(Ai - K2 - m - A) 
r(Ki - A2 - m - A)r(A2 + K2 + A - l)x 
r(Ai - (m + 1) - A)r(Ki - (m + 1) - A)r(A2 + A)r(K2 + A) 



E 



r(Ai)r(A2)r(Ki)r(K2)r(a2)r(Ai + A2 - i)r(/^i + ^2 - 1) 

(-4)-*i-*2-*3(_^)^^_^^^_^^^(_A - m)s^+s2+sz{l - A - m - (72 + si + S2 + sa)^ 



si!s2!s3! 

X (Ai - K2 - m - A)^_si_s3 (ki - A2 - m - A)^_si_s2 (Ai - (m+ 1) - A)53 (ki - (m+ 1) - /S)s2 
(A2 + A)^_S3(«;2 + A)^_S2((72 - (m + 1))52+S3(3 - A - (72 - A2 - K2)m-S2-S3 (4-23) 
where 

{x)k = x{x + I) . . .{x + k - I) 

and A is defined by Eq (4.18b) with ai = —m. 

It is interesting to consider Eq (4.23) in the hmit (72 — > —n (corresponding quantum 
field win have parameter a = —mbuji — nbuj2 and will be completely degenerate). If 
we substitute (72 = —n in Eq (4.23), then for general parameters A^ and function 
3(Ai, A2; Ki,K2; —m, —n) will be zero due to the factor l/r((72). This represents the fact, 
that the quantum three-point correlation function C{—mhu\ — n6a;2, cki, 0:2) (where = 
hr}k) with completely degenerate field and two arbitrary fields equals to zero. However, if 
one tunes the parameter ct2 for the fixed value of the parameter ai in a special way (there 
are only finite number of such possibilities), then this correlation function will be infinite, 
namely, it will have a double pole in this limit. In this case, as was explained in section 
1, it is reasonable to study the structure constants of OPE, which are defined as the main 
residues of three-point correlation function. At semiclassical level it means, that one should 
fix parameters Ai and A2 and for given value of the parameter ai = —m find such values 
of the parameters ki and K2, that function (4.23) has a double pole in the limit (72 — ^ — n. 
Namely, if parameters k\ and K2 approach to the values 

= A2 - n 2/ - A;; K2 = \i - m + 2k - I, (4.24) 

i. e. r/2 — ^ ?7i — nui — mu2 + lei + ^^2) then the double pole appears for integer k and I 
restricted by the conditions 

A;>0; l<m + n; I > k - m. (4.25) 
Prom the quantum point of view, it means, that correlation function 

C{—mbu)i — (nb + e)L02, ai,ai — nbcui — mbuj2 + Ibei + kbe2) ~ ^ (4.26) 

becomes infinite in the limit e — > 0. It follows from the fact, that this correlation function 
coinsides up to Weyl transformation with correlation function 

C{-mbLOi - {nb + e)u2, ai,al- bh^J = 

= R-\ai - bh^JC{-mbLOi - {nb + e)uj2, ai,2Q-ai+ bh^J, (4.27) 
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where 

^mn = muJi + nu)2 — kei — le2, hmn — 'muj2 + nuji — ke2 — lei (4.28) 

and R{a) is the maximal reflection amplitude defined by Eq (1.29). Correlation function 
in the r. h. s. of Eq (4.27) has a double pole, because the sum of all parameters satisfies 
the screening condition in the limit e — ^ 

-mbuJi - nbuj2 + ai + 2Q - ai + bh^^ = 2Q - kbei - lbe2- (4.29) 

The main residue in this pole should be associated with structure constant 

C{-mbui - {nb + e)u2,ai,2Q - ai + bh^J,J = ^0^^!^";^^,^^^^,. (4.30) 

So, in this case it is reasonable to consider the semiclassical limit of this structure constant. 
Maximal reflection amplitude R{bri), which due to Eq (1.28) coinsides with inverse two- 
point correlation function, has the following semiclassical limit (here rj = kiuji + K2UJ2) 

(,2\ 2ki+2k2 
^ («i-1)(k2-1)(«i + «2-2). (4.31) 

Semiclassical limit of the structure constant (4.30) can be obtained by multiplying Eqs 
(4.23) and (4.31), substituting (4.24), taking the limit o"2 —n and finding the main 
residue in this limit (one should take also into account the factor (7i-^5^)~-^i~-'^2-«i-«2-o-i-ct2 
in Eq (4.14), which relates function (4.21) with the semiclassical limit of the three-point 
correlation function). We see from Eqs (4.14), (4.27) and (4.31), that the structure constant 
defined by Eq (4.30) has a smooth semiclassical limit independent on partition function 
Zq. If we parameterize a = bXiuJi + b\2UJ2, then the semiclassical limit of the structure 
constant C'"^6t;T-n6a)2 a defined by the relation 

CTr^LT-nU.,,a^,^ln{>^^M (4-32) 

can be written in a form 



Cmn(-^l)'^2) 



TT nb'^ \'^~^'' mini 



r(Ai)r(A2)r(Ai + A2-i) 



Sm'n(Al, A2) 



(4.33) 



r(Ai -m + 2k-l- l)r(A2 -n + 2l-k- l)r(Ai + X2-m-n + k + l-2) 
with function S^^„(Ai,A2) defined as 

xT{Xi + \2-m + k-l)r{Xi-m + k-l-l)r{X2-m-n + l-l)r{Xi + X2-m-n + l-2)x 

r ('_4')-'Sl-'S2-S3 

X } , , , {l-l+k-m+n+si+S2+S3)siik-l-m)si+s2+s3{-k)m-si-s3'X 

^ A_ Si!s2!s3! 

•51,'S2,S3=0 

x(-m)si+s2+s3(^-"^-'^)m-si-s2 (-l-"^-'^)s2+s3 (Ai + fc-m-l - 1)^3 (A2-F^-m-n- 1)^2 X 



X (Ai - m -I- k)jn-si (A2 + Z - k)jn-sz (3-A;-|-m-|-n-Ai- A2)m-S2-S3 



. (4.34) 



-36- 




7 10 k 



Figure 2: Fusion rules for the completely degenerate field specified by the highest weight = 
mu)i + nu)2 with {m,n) = (7,3). Each point corresponds to admissible pair {k,l) in Eq (4.33). 
Points on the boundary correspond to the weights with multiplicity one, points on the next layer 
boundary (which arc shown by circles) correspond to the weights with multiplicity two, points, 
shown by crosses, correspond to multiplicity three and points, shown by stars (hexagon is degenerate 
into triangle), correspond to multiplicity four. Total number of points (including multiplicities) 
coincides with dimension of the representation (7, 3): (m + l)(n + l)(m + n + 2)/2 = 192. 



From the explicit form of the function Ilm,n{^i, ^2) it can be shown, that it is non-zero 
only if numbers k and besides inequalities (4.25), are restricted by the conditions'^ 

/ > 0, k<m + n, k>l-n. (4.35) 

It is convenient to picture admissible pairs {k, I) (they are defined by the conditions (4.25) 
and (4.35)), as a set of points on a plane. Namely, function (4.33) is non-zero only if points 
{k,l) lay inside of a hexagon 

k>0; l>0; m + n>k; m + n>l; n + k>l- m + l>k. (4.36) 

An example of fusion rules for the completely degenerate field specified by the highest 
weight O = muji + nuj2 with (m, n) = (7, 3) is shown on a fig. 2. We note also, that these 
fusion rules coinside exactly with quantum fusion rules, which were defined in section 1. 

The sum in (4.34) can be reduced to a simple product for all points (A;, Z), which have 
multiplicity one (i. e. for points on the boundary of the hexagon). It follows from the fact, 

These conditions follow also evidently from the functional relations (4.37) (see below). 
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that function E^„(Ai,A2) satisfies remarkable functional relations^"^ . Namely, 

^m!n(-^li-^2) = 5]^*j„(A2, Ai) = 

= + k-mM = ^i%i-k,n+k-li>^i -l + k,X2-k + l). (4.37) 

In quantum case the structure constant, which corresponds to the weight with multi- 
plicity 1, as was noticed in section 1, can be represented as a product of 7— functions. 

Using functional relations (4.37) it can be shown, that the number of terms in the sum 
(4.34) can be reduced to the minimum of numbers 

N{m), N{n), N{k), N{1), N{m+n-k), N{m+n-l), N{m+l-k), N{n+k-l), (4.38) 

where N{m) is given by Eq (4.22). It is evident from the figure 2 that this number 
depends only on the multiplicity [} of the corresponding weight of the representation 
with highest weight mui + nL02. In quantum case expression for the structure constant 
(4.27) for completely degenerate field can be reduced to the 4({) — l)-dimensional Coulomb 
integral [22]. For f) = 1 this correlation function can be given, as a product of 7- functions. 
For () = 2 it can be expressed in terms of hyper geometric function of the type (3, 2) at unity, 
while for f) > 2 it has more complicated structure. This fact clarifies the statement, which 
was done in section 1, that the complexity of the structure constant depends drastically 
on the multiplicity of the corresponding weight. 

We note also, that due to Eq (1.49) structure constant (4.30) coinsides up to multi- 
plicative factor with Coulomb integral 

CTr^Sr-nh^^^oc = {-T^lJif^^hA-rnhoJi-nhu2,a,2Q-a+mhui+nbu2-khex-lhe2), (4.39) 

where integral Ik^i{—mhijJi — nbi02, 2(5 — a + mbuji + nbi02 — kbei — lbe2) is defined by Eq 
(1.33) for the case n = 3, which is studied in details in appendix C. Using the notations of 
the appendix C this integral equals 

Ik^l{—mbu\ — nbu)2-, a, 2Q — a + mbui + nbi02 — kbe\ — lbe2) = 

= ^^feK-^''Al,-62A2,-m6^-n62), (4.40) 

where the integral Tki{—h'^\i, — 6^A2, —mb'^, —nb"^) is given by Eq (C.l). Using asymptotic 
(C.5) for this integral, we can obtain additional representation for the semiclassical limit 
of the structure constant (4.30). 

Let us say few words about the semiclassical limit of general 5i{n) TFT. This limit is 
governed by the classical equation^^ 

n-l 

ddcf) + J] efce('='='*) = 0. (4.41) 
k=l 



^''We suppose to give a proof of quantum version of relations (4.37) in Ref [22]. 
^^In Eq (4.41) we have set for shortness 7r/i6^ = 1. 
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One can show, that as consequence of Eq (4.41), field e (^i'*^) satisfies holomorphic and 
antiholomorphic differential equations of the order n 



(-5" + T9"-2 + . . . )e-(^i'<>) = 
(-5" + Jd''-'^ + ... )e-^^^^'^^ = 



(4.42) 



where by . . . denoted terms with derivatives of the smaller order. As a consequence of 
Eq (4.42) field e"*^ can be presented as a bilinear combination of the holomorphic and 
antiholomorphic solutions to (4.42) 



fc=l 



It follows from Eq (4.41) that fields e ^^'',<t') f^nh. k ^ 1 have a form 



(4.43) 



g-(a;fc,</.) ^ 

h<32<-<jk 



- ]1 "^32 



(9* 



3k 



31 



"^32 



9^31 9^3 



32 



32 



3k 



3k 



(4.44) 

The consistency condition is that the product of Wronskians of holomorphic and antiholo- 
morphic solutions equals to unity 



*n] =1. 



(4.45) 



In the case of light exponentials, all currents in Eq (4.42) are equal to zero. In this case 
functions will be the polynomial of degree n — 1 



1^ 1^ 2 n - 1 



(4.46) 



The condition (4.45) transforms to SL{n,C) constraint for the matrix aj^ 



det 



(")\ 
1 * 

in) 
2 



\ (1) (2) (n) I 

\Cln Oifi ■ ■ ■ dn I 



(4.47) 



Semiclassical limit of the correlation function of the light exponentials described by the 
integral 



1 



(Hr?i {zi,Zi) . . . Vbr^j, {ZN, Zn)) 



N 

n 

k=l 



Xvk,Hzk,Zk)) 



(4.48) 



here dil{a~^'^) is SL{n,C) invariant measure. We suppose to consider semiclassical calcula- 
tions, which were done here, for n > 3 in other publication. 
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Figure 3: Cylinder used for minisuperspace approximation. 



5. Minisuperspace limit 

It is interesting to consider another limit of sl{n) TFT at 6 — > 0. In this limit in the 
Hamiltonian picture, associated with radial quantization^^, we take into account only the 
zero mode dynamics (minisuperspace approach). In this approximation the state created 
by the operator Vg+jp^. corresponds to the wave function 

where x is a zero mode of the field ip. The function '^^p\x) {si{n) Whittaker function) 
satisfies Scrodinger equation 

(^-Vl + 27r//|^e^('^'")j ¥^'\x) = PH^^\x), (5.1) 

and in the region (cj, z) < (Weyl chamber) possesses the asymptotic 

^>f{x) ~ exp(i(P,x)) + ^ 5,(P)exp(i(s(P),x)), (5.2) 

where the sum runs over all elements of the Weyl group W besides identical. The coefficients 
Ss{P) are known exactly [38] and can be obtained from the reflection amplitude (1-17) in 
semiclassical limit 6^0 



(- 



i{s{P),e) 



e>0 



One can show, that conditions (5.1) and (5.2) determine the Whittaker function unambigu- 
ously. The minisuperspace approximation is valid if Pj/h are fixed at the limit 6^0. If we 



^®Ono should take the geometry of semi-infinite cyUnder of circumference a € [0, 27r] (fig. 3) and consider 

the states on the circle. 
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take 03 = ibq and P/. = bpk for k = 1,2, then the minisuperspace Umit of the three-point 
correlation function (1-31) can be represented by the integral 



C{Q + ibpi,Q + ibp2,ibq) J *g(x)M/g(x)e^''(«'^) . (5.3) 

The theory of the s[(n) Whittaker functions has some long history [39]. In particular, 
different explicit integral representations for these functions exist [40, 41] (wc will use 
here representation, which can be extracted from the paper [41]). All these functions 
are build from the Macdonald function Ku[y) (which can be given by the integral (D.l)) 
by the recursive integral representation. To make sense the future statements we define 
*p^(a;) = ^'p^(a;) = 1. It is useful also to introduce the variables 

Vk = b-'^e'('>"-'>/^, (5.4) 

and new function ^p^(yi, . . . , yn-i) through the relation 

e>0 '^-^ 

(5.5) 

The recursive relation connects function ^'p ^ with function *p, ^-^ 



X 

k=l '- 



h tn-2 



(5.6) 



where by definition to = 0. Vector P' in Eq (5.6) defined in a following way. If vector P has 
a components Pi, P2, • • • , Pn-i in the basis of fundamental weights of the Lie algebra sl{n) 
(i. e. (P, Cfc) = Pfc, where are the simple roots of 6l(n)), then vector P' has components 
P2, P3, . . . , P„_2 in the basis of fundamental weights of the Lie algebra s[(n — 2) (i. e. 
(P',efc) = Pfe+i, where are the simple roots of sl{n — 2)). To clarify Eq (5.6) and 
definition of the vector P' we give an expression for the function *p''(yi, 2/2, 2/3) in the 
appendix E (Eq (E.l)). 

It is useful to introduce also Whittaker function in the momentum representation 

¥^\q) = J ^>^^\x)e'^i'^Ux. (5.7) 

The integral (5.3) describing the asymptotic of three-point correlation function transforms 
to 

C{Q + ibp,,Q + ibp2,ibq) ^ J dq'¥^^{q')¥^l{q-q'). (5.8) 
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Let us consider the simplest examples of sl{2) and sl{3) TFT here. For the s[(2) case 
(Liouville theory) we derive from Eq (5.6) 



Fourier-transformed Whittaker function can be easily find using Eq (D.3) 

*£)(M = -!-(^) ^ -A^ (5.9b) 

bV2^b^^ r {-ipV2) ^ ' 

The integral (5.3) in this case can be evaluated using the formula (D.4) 
I dx^'g;(x)*g)(ar)e^''^/^- = 



. (5.10) 



- r(-...!^)r(-.P.V^) (^) " , 

This function coincides with minisuperspace limit of the three-point correlation function 
for the sl(2) TFT. 

In the 5t(3) case we obtain from Eq (5.6) the expression for Whittaker function^^ 

X yi'^'^^-^^)ifi(p,eo)(2yi0TW)^i(p,eo)(2y2v/r^ (5.11a) 

Fourier-transformed Whittaker function can be easily found using Eqs (D.3) and (D.5). 
The result is expressed again in terms of gamma-functions [38] 

^(3)^. ) = 1 (^)--(^+^-^^ nLi r(z(9,a^i) + z(p, h,))T{t{q,U2) - z(p, h^)) 

^^^^^ 62V3ne>orH(p,e)) r(%p)) 

(5.11b) 

The integral (5.3) is much more complicated in this case. It is better to write it in the 
momentum representation. As a result, for the asymptotic of the three-point correlation 
function in s[(3) TFT we obtain Barnes-like integral 

C{Q + ibpi,Q + ibp2,ibq) ^ ^"^^ 



xT{i{q-q',uJi)+i{p2,hk))T{i{q-q',uj2)-i{p2,hk)) ■ (5.12) 



'^This function was firstly derived in Ref [42]. 
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This integral can be calculated exactly in terms of gamma functions if g = scoi ot q = suj2, 
as was first noticed in [43]. 

In the case of higher n, Whittaker function is more involved object. The problem of 
finding the Fourier transform of the product of two Whittaker functions was considered in 
[44]. The most general situation, when the answer can be expressed in terms of Gamma 
functions, is g = suji oi q = .sc(j„_i. The generalization of the explicit sl{2) formula (5.10) 
for the case of general n has a form 



We note that expression (5.13) coincides exactly with the minisuperspace limit of the 
three-point function (1.39). 

Conclusion 

In this paper wc have considered in details particular examples of three-point correlation 
functions (Kn (^ii ^i)Kt2 (^2) ^2)Kk3(^3) ^3)) in s[(n) conformal Toda field theory, which can 
be expressed in terms of known in mathematics special functions. If any vector parameter 
ai, a2 or 0:3 is proportional to the first or to the last fundamental weights {ui or Un-i) of 
the Lie algebra s[(n), for example 03 = then three-point correlation function can be 

expressed in terms product of so called T-functions (see Eq (1.39)). Unfortunately, general 
situation is much more complicated. For example, if one shifts slightly parameter 
>c(jOn-\ — hu)i then corresponding three-point correlation function can be expressed only in 
terms of Coulomb integral (see Eq (2.47)) or equivalently in terms of higher hypergeometric 
functions (2.48). It is difficult to expect that general three-point correlation function can 
be expressed in terms of known functions. 

As we see from the results of sections 3, 4 and 5, where the semiclassical and minisu- 
perspace approaches to TFT were considered, three-point correlation function is already 
nontrivial in these cases. For example, in "heavy" semiclassical limit, developed in section 
3, a problem of finding it is rather difficult due to the presence of accessory parameters. 
These accessory parameters disappear only for the case of the Lie algebra sl(2), which 
corresponds to the Liouville field theory. This is one of the reasons why the three-point 
correlation function can be found in quantum LET exactly. For sl(n) TFT with n > 2 
there is no (to our knowledge) simple regular procedure to obtain accessory parameters. 
In the "light" semiclassical limit (section 4) and in the minisuperspace limit (section 5) 
it is possible to derive the expression for the three-point correlation function in terms of 
finite dimensional integrals. Generally speaking it is not evident, that in quantum case it 
is also true. The only thing, which can be done, is to find all cases when the quantum 
three-point correlation function can be expressed in terms of finite dimensional integrals. 
We will study this problem in the second part of this paper [22] . 




(5.13) 



e>0 
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A. The Coulomb integrals 

Here we will discuss the problem of calculation of the 5l{n) Coulomb integrals. They appear 
in the theory of massless n — 1 component scalar field (p, as expressions for the correlation 
functions of the exponential fields Va = e^'^''^\ We will concentrate ourself on three-point 
correlation functions 

n-l 

/.i....„_i(«i,a2,a3) = (F„,(oo)K«.(l)^a3(0) J] (^-1) 

k=l 

here Qk is a screening field Qk = J e''^'^'^''^^(Pz and are some non-negative integers. 
Correlation function (A.l) is non-zero only if the screening condition 

n-l 

ai+a2 + a3 + b'^ s^efe = 2Q 
k=l 

is satisfied. In this case correlation function (A.l) can be rewritten using the Wick rules 
^i...s„_i(ai,a2,Q;3) = 

/n—1 n—2 
fc=i /=i 

where Vs^{tk) is defined by 

T^sSk) = f[\tf-t'i\\ (A.3) 

In Eq (A. 2) we have used the notations tk = (t^^\ . . . ,ifc*'°^) with t^^^ being the coordinate 
of the j-th screening eK<^k,v) and we denote 

(A.4) 

\tk-k\ = f{f{\tf-t\^\ iik^i 

1=1 j=i 
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We will study particular case of the integral (A. 2), which corresponds to the value of 
parameter = >CLOn-i 

«^i...s„_i(ai • • • an-i|c) = 

/n—l n—2 
k=l 1=1 

with 

c = -6x, Ofe = -6(a2, ejfc). (A.6) 

This integral can be calculated using the following identity between integrals of the dimen- 
sion 2m and 2n [45, 29, 30] 



n n+m+1 



dHn{w)Vn{w)\[ W 

1 1 



i=i j=i 

n+m+l 

n +Pi) . ni n+m+1 

j 2<J J— .7 — 



I _ |-2pj-2 



(A.7) 

where T>n{t) is defined similar to Eq (A. 3) and equals 

^n{t) = '[[\ti-tj\\ (A.8) 

i<j 

Measure of integration is defined similar to Eq (A. 4) and equal djJLniw) = YYj=i d'^Wj. 

Below we list the main steps of calculation. Using integral relation (A.7) function (A. 5) 
can be calculated as follows 

• First, it is convenient to represent the factor V^^^^ {ti) in Eq (A. 5) as 

and substitute the factor 'Dj^~'^^^{t]) using Eq (A.7) with n = si — 1 and m = 

-^^^-''^ti) = / P.,_i(yi) \yi - h\-'''-'d^^s,.,{y,). (A.9) 

Note that the number of variables y\ is equal to si — 1. 

• Second, the integral over variables ti should be converted using Eq (A.7) to the form 

_ 7^l-"^-^(-b^)7(l + ai) |-2fe2+2ai ir . \2+2b^+2a^ j.l+2b^ (. ^ ^ 

-7(2 + (.2-.i + l)62 + ai) 

X v-'s^"^{y^) J P,,_i(y2) Im - ir'"''' \yi - m?'" Im - h\-'''-'df,s,-i{y2). 

(A.IO) 
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The number of integrations over variables y2 is equal to S2 — 1. We note that factor 
V^^Zi^ (yi) in the r. h. s. of Eq (A. 10) combines with a factor in Eq (A. 9) 

to the standart combination T)~g^^^{yi) and interaction between si — 1 points yi and 
S2 — 1 points y2 has a standart form \yi — y2 



Third, we note that the factor Vl^"^'^^ (t2) , appearing after the second step, combines 
with the factor V~^^^ {12) in Eq (A. 5). Hence we can take the integral over variables 
t2 in a way similar to the second step 



7(3 + (S3 - S2 + 2)6^ + oi + 02) ^ 

(A.ll) 

• Repeating this procedure we will lower the number of integrations at every step. The 
last integral over variables tn-i will be different from the integrals appearing at the 
previous steps. Namely, 

/ P.„_l(U-l) |4-lP^ |4_,-l|2(-2)(m^)+2Ea. |C-yVl|-^'''-^d/X.„_,(t„_l) 

= ^^n-i-i(_^2N 7(1 + chin - 1 + (n - 2)6^ + ai + • • • + an-i) 
7(n + c + ai-\ h a„_i + (n - 1 - Sn_i)62) 

X 25r;_f ''^2/n-l)|yn-lP^-^'''bVl - l|2H-2{n-3)(l+;.^)H-2K+...+a„_0 (^.12) 

• As result, we obtain the recurrent relation 

Jsi,...,s„-i{ai,a2, . . . , a„_i|c) = 

= K{ai,a2, ... , an-i\c) Jsi-i,...,s„_i-i(ai - 6^,02 . . . an-i|c - 5^) (A. 13) 

with 

7(-si62) 7(1 + c)7(n - 1 + ai + ■ ■ ■ + On-i + (n - 2)6^) 



K{ai,a2, . . .,an-i\c) 



7"-i(-62) j(n + c + ai + --- + a„_i + (n - 1 - Sn_i)62) 

,,ff 70- + «i + --- + % + (i-i)^'^) 

7(1 + i + «1 + • • • + + - Sj + i)62) 



We note, that if we substitute parameters and c from Eq (A. 6), then the solution to the 
recurrent relation (A. 13) gives the expression for the integral (1.38). We note also, that 
recurrent relation (A. 13) can be used to continue integral Jsi,...,s„_i (ai, 02, • • • jOn-ijc) to 
the non-integer values (it gives the expression for the three-point correlation (1.39)). 
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B. Simplification of the integral (4.18) 
We start with the integral (4.18) 



a(Ai,A.;.i,K.;ai,a.) =4-+-y ^a^a. -j-d-^adHd^ (B.l) 

where is given by Eq (4.18a). First, we use Feinman representation 

1 _ r(t;) f ^Ai^A2^Ki drd^dt 

z^'z^^z^'z^' ~ r(Ai)r(A2)r(Aci)r(K2) J {Z4 + z^t + + z^r)^ VJJ 

with 

f = Al + A2 + + K2- 

After that, we can calculate integral over d'^a, d'^b and over dg with the result 

r(Ai)r(A2)r(«i)r(«:2) J ^ ^ ^ * % r ^ t ^^-^^ 

with 

r = l + t + T + ^ + T|cp + t|c + z^p, f = (1 + + iri/^; 

P = I(l + t(l + |c + .n) + |(C + r(l + |cn)+^(l + |c+^|^). 

The problem is that the quantity P is not quadratic in the variable c. In order to proceed 
simplification, we use the following trick. We can multiply our expression (B.3) by 

1= / _ / dps'P= ds5{s-l) (B.4) 
J— 00 ^ J— 00 J— 00 

and insert s somewhere into (B.3). More exactly, we will need to do that four times. We 
just show the places of insertion of different Sfc 

r ^ .S4 + t{1 + |cp)) + (1 + i(l + |c + z^P)) , f ^ (1 + + sstri/^; 

The integrals over Sk can be calculated exactly (first over si, second over S2 etc). After 
that, the integrals over t, r and ^ will be of the type (D.5) and also can be calculated. The 
remaining integral over dv and d^c will be 

y'^2(5+A-.i-.2)(l^|^|2)3i+.2-.3-Ai^^^|^^^|2^S2+.3-.H^^^|^^^|2yi-5 d^^2^ 

Using technique, described above, one can reduce integral (B.5) to one dimensional integral. 
As result, the integral (B.l) can be reduced to five dimensional Barnes-like integral. By 
using the first and the second Barnes lemmas (D.6) and (D.7), one can reduce it to three 
dimensional integral (4.19). 
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The integral (4.19) can be also rewritten in a different form in terms of Tricomi func- 
tions, which are defined by the integral representation 



■^{a, c\x) = -— / dt e-^*r-^ (1 + t) 

Jo 



c—a—l 



(B.6) 



This function can be expressed through the confluent hypergeometric function 

4(a,c|x) = l+-x + 5^A_|x= + ... (B.7) 

as 

c\x) = Ha, c\x) + llxi-'^ $(a _ c + 1, 2 - c\x) (B.8) 
L [a — c + 1) \. [a) 

with *(a,c;0) = y^^-^^i) • Tricomi function satisfies the following relation 

^-(0, c|x) = x^-" *(a - c + 1, 2 - c|a;). (B.9) 
The integral (4.19) can be rewritten as 

5(Ai, A2; Ki,K2;ai,a2) = 4^i+«i+-i-^x 

r(Ai + Ki + (Ti - A - 2)r(A2 + K2 + as + A - 2) 



r(Ai)r(A2)r(Ai + A2 - i)r(M)r(K2)r(Ki + k2 - i)r((7i)r(c72)r(ai + ^2 - 1) 

X J dudsdy ^Ai-l^i _^)Ki-l^Ai+Aci-A-2^^ _^)A2+«2+A-2g-5^Ai+Ki-cT2-A-l^ 

X Fi(4t(l - t)s) F2{sut) F^{su{l - t)), (B.IO) 

where 

Fi{x) = r(c72)r((72 + A) *(a2 + A, 1 + A|a;), 

F2{x) = r(Ki +K2- l)r(ai + Ki - A - 1) *(ki + K2 - 1, 1 - ai + K2 + A|x), (B.ll) 

F3(x) = r(Ai + A2 - l)r(ai + Ai - A - 1) ^'(Ai + A2 - 1, 1 - ai + A2 + A\x). 

This form of the integral (4.18) is very convenient to obtain its limit at ai — —m and 
(72 —n considered in section 4. 

C. Properties of the sl{3) Coulomb integral 

In this appendix we study the properties of the sl(3) integral 

i=l j=l i=l J=l 

(C.l) 
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where T)k{t) is defined by Eq (A. 3). Using the integral identity (A. 7) one can show, that 
function Ik,l{c(i,ce2, (^i, P2) satisfies the set of functional relations, which are generated by 
two basic relations (we suppose, that I > k): 

Ik,l{(^i,a2,Pi,P2) = S^^](Q:i,Q:2,/3i,/32)2rfe,K«i>/3i,a2,/32), (C.2) 

where Pi = Pi + {I — k)b'^, 0.2 = 0.2 — {I — k)b'^ and 

l—k—l , .,2 



i-k-i , .^2^ 



7(2 + ai + a2 - (i - 1)62) ^ 7(2 + /3i + /32 - (i - !)«'') 



X 



-j-j- 7VZ, -r m -r «2 - v:/ - -1-;^^ ; j-j- 



7(2 + ai + A - (i - 1)62) 7(2 + 52 + ^2 - U - 1)62) 
and by the relation 

f2) 

^k,l{(^i,o:2,Pi,P2) = 'Blj{ai,a2,Pi,P2)x 

X Ik,i{ai, -2-ai-a2 + {l- 2)6^, (3^^-2 - p^ - P2 + {I - 2)6^) (C.3) 

with 

"m("1'"2,/?1,/?2) = 

^ 7(1 + ai - i62)7(l + ^1 - j62)7(2 + cti + Q2 - (j - 1)62)7(2 + P1 + P2- {j - 1)6^) 
J=o 7(2 + ai + /3i-(i-A;-l+i)62)7(3 + ai + a2 + /3i + /32-(A;-2 + i)62) 

Relations (C.2) and (C.3) can be used for the analytical continuation and sometimes for 
the simplification of the integral (C.l). 

Integral (C.l) can be calculated exactly if A; = or Z = and also if one of the 
parameters ct^ or Pk equals to zero (sec Appendix A). In the case k = 1 (or Z = 1) it also 
can be reduced to known functions. To show it we apply integral relation [29, 30] 



1-2 
=n 



7(-(j + 2)62) 7(1 + a2- i62)7(l + P2 - jb^ 



.^0 7(-62) 7(2 + a2 + /32-(^-l + j)62) 

X 1 y \u\2a2-2{l-l)b^ 1^ _ ^^2/3,-2{l-l)b^ |^ _ t\2lb^^2^_ (c.4) 

Relation (C.4) allows to reduce integral (C.l) to the four-dimensional integral 

|i|2"l \t - l|2/?l \u\2a.-2{l-l)b^ 1^ _ -L|2^,-2(;-l)b2 |^ _ ^^2lb^ ^2^ ^2^^ 



which can be expressed in terms hypergcometric function of the type (3, 2) using Eq (2.48). 
For A: > 1 integral (C.l) can be reduced to 4A:-dimensional Coulomb integral. We will give 
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the explicit expression for this integral in Ref [22]. Here we give two different asymptotics 
at 6 —> of the meromorphic function defined by the integral (C.l). First asymptotic is 
(we assume that I > k) 

6— ♦O 

^ (-l)HA2);-fc(AC2);-fc ^ 

{Xl + \2 + Kl + K2 + 1 - 2)fe(Ai +Ki+k-l- l)fe(A2 + K2 + l-k-l)i 



k 



X X /l-'Sl-S2-S3 



{-k)si+S2+S3{-l)si+S2+S3i^ - 2fc - Aj - + gj + §2 + ^3)51 



y 4- 

Sl,S2,S3>0 



X {>^l)k-si-S3{l^l)k-si-S2{^l +Ki + k-l- 1)52+53(3 - Ai - A2 - Kl - K2 - k)k-S2-S3X 
X (ki +K2- l)s2('^l + A2 - l)s3(^ -k + A2)fe-S2(^ -k + K2)k-S3- (^-5) 

Second asymptotic is 
Ik,i{-1 - Al6^ -1 - A26^ -1 - Kib^ -1 - AC26') ^ (-p)*^^' X 

3^0 s^O ' (Ai)fc_,i(A2)z-,2(Ai + A2 - l)k-s^ {ki)s^{k2)s2{ki +K2- l)s2 ' ^ ' ^ 

where Cj, are the binomial coefficients. 
D. Useful formulae 

Here we collect some basic facts concerning Macdonald function K^{y) 

• Integral representation 

My) = -J^ jt''eM-y{t + yt)/2). (D.l) 

• Asymptotic formula 

K42y)^^{r{-,y)y^ + Ti,y)y-'^) at y^O (D.2) 

• Mellin transformation of single Macdonald function 

f.™f.J,r(^i±i^)r(^) (D.3) 

• Mellin transformation of the product of two Macdonald functions 

/■°° . dy c" ^ /A + u + z/\ /A + u-i/ 



A — /i + i^\ /A — ;U — / ^-M+i^ ^ c 



2 / V 2 



X r r ' r ^ F 



A 



1-4)- (D.4) 



a2 



here F denotes the hypergeometric function of the type (2, 1] 
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Beta-like integral 



Barnes first lemma 



dt A, 0,0 ir(4)r(-5-4) , , 



27rz J_i^ T{a + /3 + J + 5) 

Barnes second lemma states that 

1 f r(ai + s)r(a2 + s)T{a3 + s)T{l - /3i - s)r(-s)c?s 



(D.7) 



27ri 7 r(/32 + s) 

^ r(ai)r(a2)r(a3)r(l + ai)r(l - /3i + a2)r(l - /3i + «3) 
r(/32 - ai)r(/32 - a2)r(/32 - as) 

provided that /9i + /32 = ai + Q!2 + "3 + 1- 

E. Example of application of the recursive relation (5.6) 



In this appendix we explain how to use recursive relation (5.6). For example, using Eq 

>(4) 
P 



(5.6) one obtains for function ^p'*(yi, 2/2, J/a) exact expression 



/* 00 /* 00 



JO 



X K ,(P,e,) (^2y2gJ i^ »(p^eo) (^2yi-^(l + tr^)J X 

X i^ .(P,en) (^22/2-^/(1+ i!)(l+i2"')) i^ »(fien) (^2^3-^/(1 + *!)) (E.l) 

In Eq (E.l) we substitute 

fy2^1 = f2y2^) . (E.2) 



t2 J b \ t2 

As we see from Eq (E.2), it is convenient to think that P' = P, but vector P' lives on a 
lattice with cutted-off ends. Symbolically it can be pictured as 



ei 62 63 6„_3 6„_2 6„_i 62 63 e„_3 e„_2 

Using function (E.l) we can reconstract function ^'p^ (2/1,2/2, 2/3, 2/4, 2/5) and so on. 
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